Trimers in the resonant 2 + 1 fermionic problem on a narrow Feshbach resonance: 
Crossover from Efimovian to Hydrogenoid spectrum 
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We study the quantum three-body free space problem of two same-spin-state fermions of mass 
m interacting with a different particle of mass M, on an infinitely narrow Feshbach resonance with 
infinite s-wave scattering length. This problem is made interesting by the existence of a tunable 
parameter, the mass ratio a = m/M. By a combination of analytical and numerical techniques, we 
obtain a detailed picture of the spectrum of three-body bound states, within each sector of fixed 
total angular momentum I. For a increasing from 0, we find that the trimer states first appear 
at the Z-dependent Efimovian threshold Oc , where the Efimov exponent s vanishes, and that the 
entire trimer spectrum (starting from the ground trimer state) is geometric for a tending to ai 
from above, with a global energy scale that has a finite and non-zero limit. For further increasing 
values of a, the least bound trimer states still form a geometric spectrum, with an energy ratio 
exp(27r/|s|) that becomes closer and closer to unity, but the most bound trimer states deviate more 
and more from that geometric spectrum and eventually form a hydrogenoid spectrum. 

PACS numbers: 34.50.-s,21.45.-v,67.85.-d 



I. INTRODUCTION 

The quantum few-body problem is the subject of a 
renewed interest [l[ , thanks to the possibility of experi- 
mentally studying this problem in a resonant regime with 
cold atoms close to a Feshbach resonance |3 • In this res- 
onant regime, the s-wave scattering length a associated 
to the interaction among the particles can be made much 
larger in absolute value than the interaction range. This 
has in particular allowed to study the Efimov effect in the 
lab, up to now for three bosons and for three distinguish- 
able particles 0^, that is the emergence for 1/a = of 
an infinite number of trimer states with an accumulation 
point at zero energy in the vicinity of which the spec- 
trum forms a geometric sequence. Whereas the existence 
of an infinite number of bound states is common for long 
range interactions (vanishing for diverging interparticle 
distance ri2 as l/r^2 ot more slowly), this is quite in- 
triguing for short range interactions. Initially predicted 
by Efimov for three bosons, this effect can actually take 
place in more general situations [J], in particular in the 
so-called 2-1-1 fermionic problem if the extra particle is 
light enough [5j. 

What we call here the 2 -fl fermionic problem consists 
in the system of two same spin state fermions of mass m 
interacting with a particle of mass M of another species. 
It is assumed that there is no direct interaction among 
the fermions, whereas there is a resonant interaction be- 
tween each fermion and the extra particle, that is with an 
infinite s-wave scattering length, 1/a = 0. Furthermore, 
it is assumed that this resonant interaction is due to an 
infinitely narrow Feshbach resonance, that is of vanishing 
van der Waals range 6 — > and finite effective range rg . 

This concept is most easily understood in a two- 
channel model. In the open channel, the particles ex- 
ist in the form of atoms, and have a weak, non-resonant 
direct van der Waals interaction corresponding to the 



background scattering length abg ~ b and the interaction 
range b. In the closed channel, the particles exist in the 
form of a bound state of a fermion with the other-species- 
atom, the so-called closed-channel molecule with radius 
w b. Due to a coupling A between the two channels, 
that we shall precisely define later, the closed-channel 
molecule is coherently converted into a pair of atoms in 
the open channel, and vice-versa. For an appropriate 
Zeeman tuning (with a magnetic field) of the bare energy 
£"1110! of the closed-channel molecule with respect to the 
dissociation limit of the open channel, the s-wave scat- 
tering length a between a fermion and the other-species- 
atom is infinite. In this case, in simple models, the ef- 
fective range rg is the sum of two contributions [6| . The 
first one is, as expected, of the order of the van der Waals 
length b. The second one is induced by the interchannel 
coupling; it is expressed as — 2i?* [7], where the Feshbach 
length R^ is positive and scales as 1/A^: 



i?* = 



2, ,2 



A^fi 



(1) 



where /i is the reduced mass of a fermion and the other- 
species particle. When the interchannel coupling A is 
very weak, this second contribution dominates over the 
first one, i?* 3> 6, and this is the narrow Feshbach res- 
onance regime. An example under current theoretical 
and experimental investigation is the case of the inter- 
species Feshbach resonances of the fermionic ®Li and the 
fermionic *°K, which are narrow: The Feshbach length 
i?* exceeds 100 nrn, whereas the van der Waals length is 
a few nanometers [8l-[l0|| . 

To obtain the infinitely narrow Feshbach resonance 
model, here for 1/a = 0, one takes the mathematical 
limit of a vanishing van der Waals range b with a fixed 
non-zero interchannel coupling A. The s-wave scattering 
amplitude between one fermion and the extra particle for 



a relative wave vector k is then [3] : 



fk 



-1 



2/C ~T" f£ iL^ 



(2) 



This implies the absence of two-body bound states, since 
fk can not have a pole for k — iq, q > 0. The model can 
however certainly supports trimer states since the 2 + 1 
fermionic problem is subjected to the Efimov effect for a 
large enough mass ratio m/M 0, Q . 

The main motivation of the present work is to study 
the spectrum of trimers for this problem in free space 
for 1/a = 0, in particular to determine analytically the 
global energy scale of the Efimovian part of the spectrum, 
related to the so-called three-body parameter. This 
global energy scale is out of reach of Efimov's zero range 
theory Q but it was determined analytically for three 
bosons for a narrow Feshbach resonance in [ll|, Il2l ■ Here 
we shall generalize this calculation to the present 2-1-1 
fermionic problem. A second motivation is to determine 
the low-lying states of each Efimov trimer series, which in 
principle are not accurately described by Efimov theory, 
and to look for possible trimer states that are not related 
to the Efimov effect and may thus appear for lower mass 
ratios. 

The paper is organized as follows. After a presentation 
of the model and the derivation of an integral equation 
a la Skorniakov-Ter-Martirosian [l^j for the three-body 
problem in momentum space in section [TTl analytical so- 
lutions of this integral equation are obtained in limiting 
cases in section IIIII In the central section IIVI of the pa- 
per, we analytically explore the physics of the trimers: 
Of particular interest are the exact results on the global 
energy scale in the Efimovian part of the spectrum, see 
subsection IIVBI and the study of an hydrogenoid part 
of the spectrum in the Born-Oppenheimer regime, see 
subsection IIVCI An efficient numerical solution of the 
integral equation is used in section |V] to explore interme- 
diate regimes not covered by the analytics. We conclude 
in section rvn 



II. 



THE MODEL AND THE GENERAL 
MOMENTUM SPACE EQUATION 



A tractable though realistic description of a Fesh- 
bach resonance is obtained with the so-called two-channel 
models, where the particles exist either in the form of 
atoms in the open channel or in the form of molecules 
in the closed channel 0, [l4h20] . We use here the same 
free space two-channel model Hamiltonian H as in [21| 
written in momentum space in second quantized form 
in terms of the fermionic annihilation operators Ck, the 
extra-particle annihilation operators ak and the closed- 
channel molecule annihilation operators b^: 



with 

Hat — 

Hnwl = 

-^at— mol — 



_fk_ 

(27r)3 
(Pk 

(fikiCPk 



ii'kCkCk + ai^kflkflk 



l + a 



-Bk ) fej^fok 



(4) 
(5) 



i/o 



/ \(2-kW ^(^i2''^°ki+k2°kiCk2+h.c.j(6) 
(Pki<Pk2<Pk:i(PkA 



90 



-X(ki2)x(k43) 



[(27r)3]4 
X (27r)=^J(ki + ka - kg - k4)4 c], Ck^OkilT) 



Whereas the Ck obey the usual free space anticommuta- 
tion relations 



{ck,4,} = (2^)35(k-k') 



(8) 



the statistical nature (fermionic or bosonic) of the extra 
particle and of the closed-channel molecule does not need 
to be specified here, since there will be at most one of 
the particles in the state vector. Simply, the Ck and cj^. 
commute with Ck and 6k- In the kinetic energy terms of 
the atoms iJat and of the closed-channel molecule ifmoi, 
we have introduced the free fermion dispersion relation 
i?k = h'^k'^ /{2m) and the mass ratio of a fermion to the 
extra particle: 



m 
M' 



(9) 



The internal energy E'moi of the closed-channel molecule 
is counted with respect to the dissociation limit of the 
open channel and is experimentally adjusted thanks to 
the Zeeman effect by tuning of the external magnetic 
field. i?at-moi represents the coherent interconversion of 
a closed-channel molecule into one fermionic atom and 
the extra particle, due to the coupling between the closed 
channel and the open channel. It involves the interchan- 
nel coupling constant A and is regularized by the mo- 
mentum space cut-off function %, assumed to be real and 
rotationally invariant, that tends to one a zero momen- 
tum, and that rapidly tends to zero at large momenta 
with a width 1/6, where the interaction range b is of the 
order of the van der Waals length. Note that the argu- 
ment of the cut-off function x is the relative wave vector 
of a fermion (of momentum ftka) with respect to the extra 
particle (of momentum fiki): 



ki2 



'k2 _ k]_ 

^^ m M 



aki 



where 



M 



mM 
m + M 



(10) 



(11) 



H = -ffat + -ffmol + ^at-mol + -^opc 



(3) 



is the reduced mass, so as to preserve Galilean invari- 
ance. Finally, -ffopon niodels the direct interaction be- 
tween atoms in the open channel, in the form of a separa- 
ble potential with bare coupling constant go and the same 



cut-ofF function x as in iJat-moi- This direct interaction 
is characterized by the so-cahed background scattering 
length Cbg. The inclusion of both TJat-moi and ilopcn 
allows to recover the usual expression for the scattering 
length as a function of the magnetic field B [2] , 



a{B) = flbg X 1 



AB 
B-Bo 



(12) 



if -E-inoi is taken to be an affine function of B. The quan- 
tity AB is the so-called magnetic width of the Feshbach 
resonance. 

We now derive from the model Hamiltonian a mo- 
mentum space integral equation a la Skorniakov-Ter- 
Martirosian 13| for the three-body problem of two 



fermions and one extra particle, closely following ref- 
erence [21[ downgraded from the 3 -|- 1 to the 2-1-1 
case. In the search for bound states, we take a nega- 
tive eigenenergy, E < 0, and we express Schrodinger's 
equation — {H — E)\'i) for a ket of zero total momen- 
tum and being the sum of general ansatz with zero or 
one closed-channel molecule: |^) = ji/'Sat) + IV'iat+imoi), 
with 



IV'Sat) 



d^ki(fk2d^k3 



1^. 



1 at+1 iiiol 



) = 



[(2^)3]3 
^(ki,k2,k3)4^4 



(2^)3<5(ki 

Lio) 



(2^)^ 



:Bik)bl^cl\0). 



ka) 

(13) 
(14) 



Thanks to the fermionic antisymmetry we can impose 
that A(ki , k2 , ks) is an antisymmetric function of k2 and 
ks. On the contrary, a fermionic atom and a closed- 
channel molecule are distinguishable objects and there is 
no exchange symmetry constraint on the function B(k). 
Projecting Schrodinger's equation on the subspace with 
three atoms, and using £' < 0, we are able to express A 
in terms of B and of an auxiliary unknown function B 
obtained by a partial contraction of A: 



B(k3) 



d^kid^k2 



[(27r)3]2 
xA(ki,k2,k3 



X(ki2)(27r)35(ki 



(15) 



More precisely. 



A(ki,k2,k3) 



A/2 



E ~ [aEk, + E^, + Ek,) 

X [x(ki2)i?(k3) - x(ki3)i?(k2)] (16) 

where the convenient unknown function is actually D 
such that 

AD(k)==AB(k)-H2goS(k). (17) 

Plugging the expression (fT6|) of A into the definition (fT5|) 
of B gives a first important equation 



2 5ru N f d'hd'k2 



(27r)35(ki+k2+k3: 



E - (aEk, + Ek, + Ek,) 
X [x(ki2)C(k3) - x(ki3)-D(k2)]. (18) 



The second important equation is obtained by pro- 
jecting Schrodinger's equation on the subspace with one 
atom and one closed-channel molecule (in which case the 
direct open-channel interaction can not contribute): 



2AB(k) = [^,.el(k)-^„,ol]S(k) 



(19) 



where we have introduced what we call the relative en- 
ergy 



E,,iik) =E-[E, 



1 



-E. 



(20) 



This is indeed the relative energy of one of the fermions 
and of the extra particle, knowing that the second 
fermion has a wave vector k, since one subtracts in ([20)1 
from the total energy E the kinetic energy £'k of the 
second fermion and the center of mass kinetic energy of 
the first-fermion-plus-extra-particle. One expresses B in 
terms of D by elimination of B between (IT7|) and (|19l) . 
One then eliminates B between the resulting equation 
and (|18l) to finally obtain a closed equation for D: 



= 



fiDiks) 



d^kid^k2 



27Th^f[E,,,{ks)] J [(2^)3]2 

^^ X(ki2)x(ki3)^(k2) 



(2^)35(ki+k2+k3) 



E - {aEk, + Ek, + Ek,)' 



(21) 



The function / is related to the two-body T matrix for 
the scattering of a fermion and of the extra particle ^1| , 



(k/|r(e + *0+)|k,) 



2T:h^ 



x(k/)x(k.)/(e + *0+) (22) 



where the relative wave vectors k^, k/ and the energy e 
are arbitrary (the T matrix is not necessarily on shell in 
that expression). 

The general expression of / is given in [21| . Here how- 
ever, we shall concentrate on the limit of an infinitely 
narrow Feshbach resonance. We thus take the zero-range 
limit b — )> 0, in which case the cut-off function x (of 
width ex 1/b) tends to unity. It is assumed that there is 
no resonant interaction in the open channel, so that the 
corresponding background scattering length Obg is 0{b) 
and also tends to zero. On the contrary, the interchannel 
coupling A is kept fixed, so as to keep a non-zero effec- 
tive range, and -Emoi is adjusted to keep a fixed value of 
the scattering length a. In this case, the function / for 
a positive and fixed energy e > (so that fib'^e/h? -^ 0) 
simply tends to the s-wave scattering amplitude fk on an 
infinitely narrow Feshbach resonance [7] : 



f{e + tO+) = fk 



with 



fk 



1 



-1 +ik + k^R^ 



(23) 



(24) 



The Feshbach length i?, is expressed in terms of the 
width AB of the Feshbach resonance as [3] : 



i?* 



2fj,ahgf^bAB 



(25) 



where the differential magnetic moment between the 
closed and open channels is /i;, — dE^oi/dB taken 



for B = Bq. In (|24|) . the relative wavenumber k = 
(2/ie)^/^//i since the energy is positive. For a negative 
energy, the function / has the same expression ([23]) if 
one uses the analytic continuation k — i{—2fj,eY^^/h in 

fkM- 

It is convenient to represent the eigenenergy E in terms 
of a wavenumber, setting 



E = - 



2^ 



(26) 



where we recall that fi = mM/{m + M) is the reduced 
mass and q > Q. Similarly the relative energy is rep- 

cl( 

leading to 



resented by a wavenumber, £^rci(k) = —h'^qf^i{k)/{2fi) 



qi-ci{k) = 



l + 2a ;,' 

(l + a)2 



1/2 



(27) 



For an infinitely narrow s-wave Feshbach resonance with 
a scattering length a, the integral equation resulting from 
Schrodinger's equation, the equivalent of the Skorniakov- 
Ter-Martirosian equation [13| for our problem, is then 



o=[- 



i(fc)4 
d^k' 



qUk)R.] D{k) 
D{k') 



27r2 



-kk' 



(28) 



In that equation, for the sake of generality, we have kept 
an arbitrary value of the scattering length a. In what 
follows, we shall restrict to the exact resonance location 
where 1/a = 0. 

We shall also take advantage of rotational invariance 
to reduce the integral equation to an unknown function 
of a single variable only, as done in |l3j. The eigenstates 
of the Hamiltonian H may be assumed of a fixed total an- 
gular momentum of angular quantum number I. Without 
loss of generality one can also assume that the angular 
momentum along the quantization axis z is zero. The 
corresponding ansatz for D is thus 



D(k) = r,°(k)/W(fc) 



(29) 



where the notation YJ™' (k) stands for the spherical har- 
monics yj™' {9, (j)) where 8 and 4> are respectively the po- 
lar and azimuthal angles of the vector k in a system of 
spherical coordinates of polar axis z. The unknown func- 
tion /(') then depends only on the modulus k of k. We 
see that this also fixes the parity of the eigenstate to 
the value (—1)'- In the integral over k' in (|28|) . we then 
perform the change of variable of unit Jacobian, 



k'=7eK 



(30) 



where TZ is the rotation in 
composition of its inverse: 



defined by the Euler de- 



Tl-^ =TZ^{-(p)ny{9)n,^{TT 



(31) 



where the notation Tli{a) stands for the rotation of an 
angle a around the axis i G {x, y, z}. This choice ensures 
that k/A: = TZez, where ez is the unit vector defining the 
z axis. In the denominator of the integrand of (1^51) . the 
scalar product k-k' is then transformed as fce^ -K, so that 
the denominator is invariant by rotation of K around z. 
In the numerator we use the transformation of spherical 
harmonics under rotation, see (8.6-2) and (8.6-1) in 23], 
and the relation [Y;"' (6*, -0)]* 



Yr{o,cf>y. 



y,"(7^K) 



An 
21 + 1 



1/2 



E 

mi = -l 



r/"'(k)Y;"'(K) 



(32) 



which in particular allows to pull out the factors F/"' (k) 
expected from rotational invariance. The integration 
over K is then conveniently performed in spherical co- 
ordinates of polar axis z. All the terms with to; 7^ in 
(|32p vanish in the integration over the azimuthal angle 
of K. From the expression of the spherical harmonics in 
terms of the Legendre polynomial of degree I j23i] , 



>1"(K) 



47r 
21 + 1 



-1/2 



Pj(« = K-ejA'), 



(33) 



we finally obtain the reduced integral equation to be 
solved in the sector of angular momentum I, for 1/a = 0: 



0= [q 



,,,{k)+ql,{k)Rjf^'\k) 



TT 



du ■ 



Pi{u)K^ 



k^+K^ 



^kKu 

(34) 



III. ANALYTICAL SOLUTIONS IN 
PARTICULAR CASES 

Several remarkable analytical techniques are now avail- 
able to solve the three-body problem in some appropriate 
limiting cases @, [ill, [l2, IIS ISg . Whereas we do not know 
how to solve ([M)) analytically in general, it is possible to 
find solutions when there is an extra symmetry available, 
that is scale invariance. The most standard regime corre- 
sponds to the limit i?* — >■ 0, in which case our model re- 
duces to the so-called zero-range or Bethe-Peierls model, 
where the interactions are included via two-body contact 
conditions on the wavefunction |26| . Since 1/a = 0, these 
contact conditions are indeed scaling invariant , which al- 
lows to fully solve the problem [J, [23] . Because the zero- 
range model is usually solved in position space, it is in- 
teresting here to briefly show the calculations in momen- 
tum space. In the regime of interest, where the Efimov 
effect takes place, the zero-range model is however not 



well defined, and an extra three-body condition has to 
be introduced to make it self-adjoint |2^], involving the 
three-body parameter. 

So the relevant case here is R^, > 0. The existence of 
such a finite length scale characterizing the interactions 
breaks the scale invariance. Fortunately, as shown in [l2| , 
if one restricts to the zero energy case E — 0, equations 
of the type ([M)) can still be solved analytically. This 
gives access to the three-body parameter, and thus to 
the characterization of the Efimov spectrum of trimcrs. 



A. Zero-range case (R, — 0) at zero energy 

For R^ — and at zero energy q = 0, the integral 
equation ([34)) is manifestly scaling invariant: If f^^\k) 

is a solution, the function /j; (k) = /^'^ (k/X) is also a 
solution, VA > 0, and we expect that the two functions 



Z^'-* and /} are proportional. We thus seek a solution in 
the form of a power-law. 



/«(fc)-fc-(^+2). 



(35) 



The form of the exponent results from the general theory, 
see section 3.3 in [29||: For a iV-body problem, with here 
TV = 3, for s to be a directgeneralization of the exponent 
So introduced by Efimov [J] , the exponent in psp should 
be — [s-f (3A^ — 5)/2]. For convergence issues, it is simpler 
to assume in explicit calculations that s = iS^ where 
S is real. By analytic continuation, the result however 
extends to real s also, as also shown by the real space 
calculation [30|. We inject the ansatz ([55]) in ([M|) with 
q = 0, -R* = 0, and we perform the change of variable 
K = ke^ to obtain an equation for s: 



= A,(s), 



(36) 



with the function 



-iSx 



, , , (l + 2a)i/2 /-i r+--dx e 

KAs) = - h/ duPiiu) — . 

(37) 
Using contour integration and the Cauchy residue for- 
mula, the integral over x may be calculated: Setting 
9 = arccos ( ^-r^u) , so that 9 £ [0, tt], we obtain 



+ 00 



dx 



-iSx 



2tt cosh X 



l + Q 



sin(s6') 
sin(s7r) sm9 



(38) 



Successive integration over u is simplified by taking 9 
rather than u as integration variable, with a Jacobian 
that simplifies with the factor sin 9 in the denominator 
of ([55)1 . Following [5| we then parameterize the mass ratio 
by an angle i^ G]0,7r/2[: 



arcsm ■ 



M' 



(39) 



which leads to the remarkable property 

(l + 2a)i/2 



COSI^ 



1 + a 



Also using arccos[a/(l -I- a)] 



Ms) 



COSJ^- 



smz^ 



r+l' 



d9P, 



u, we obtain 



9\ sin[s9) 



svcLV J sin(s7r) 



(40) 



(41) 



This can be further simplified taking advantage of the 
parity of the Legendre polynomial, Pi{—u) = {—iyPi{u) 
for all u, by shifting the integration variable 9 by 7r/2. 
Depending on the even or odd parity of the angular mo- 
mentum / this reduces to 



Hs) 


I even , t 


sini^^ 


Ms) 


(odd 1 





smi^ 



d9P 



d9P, 



smW 
siniy 
sin 9 
siniy 



cos{s9) 
cos(s7r/2) 

sin(s6') 
sin(s7r/2) 



(42) 
.(43) 



Eq. (j43|) will be quite useful to obtain analytical results 
on the Efimovian trimer spectrum in the large a limit, 
see subsections IIV Al and IIV 51 

An explicit expression of A/(s) as a sum of a finite 
number of simple functions of s may be obtained by rep- 
resenting the function Pi (^^) as a Fourier sum, that is 
a sum of cosn9, < n < I, n even, for an even Z, and 
a sum of smn9, 1 < n < I, n odd, for an odd I. The 
coefficients are polynomials of 1/sinz^ that arc simple to 
calculate analytically from the known coefficients of the 
Legendre polynomials Pi{u) [31|. The resulting integrals 
over 9, e.g. of sin(ri6') sin(s0), are then straightforward to 
evaluate. We finally obtain the explicit formula valid for 
arbitrary parity of I: 



1 



Ai(s)=cosz/H ,, ,x , -, 

^ ' sini/cos[(s-j-/)7r/2] 



y^^ / sin[(s + n);^] ^ _^ y sin[(s - n)z/] \ 
^^ 1 s + n s — n J 



(44) 



The coefficients c„ are zero for I — n odd. For I — n even, 



1 \ (-l)('+")/2 

2 '^"■"j (Asiniyy 



(l-n)/2 



{~4:Sin^iy)''{2l-2ky. 
^ fc!(?-A:)!(^-fc)!(^-fc)! 



E 



(45) 



k=0 

where dij is the usual Kronecker delta. 

The form (|44)l has the interesting feature that the co- 
efficients are s-independent, which makes the numerical 
evaluation of Ai(s) as a function of s particularly effi- 
cient. As a test, it is however interesting to compare to 
the transcendental equation for s obtained by the direct 
calculation a la Efimov in position space. As detailed in 
the Appendix \^ introducing the hypergeometric func- 
tion 2^^! as in ^32j to solve some differential equation, we 



generalize the formulas of [32| to an arbitrary mass ratio, 
a generalization that was done already with the adiabatic 
hyperspherical method in |33j : 

A,(s) = cos ^y + (-!)' sin' z/- ^ 2 ^ ^ 2 ^ 



X2i^l 



1 + 1 + s l + l 



* , 3.0 

,i-\ — ; sm v 

2 '2' 



(46) 



where T is the Gamma function. A variant of Eq. ((46]) 
will be quite useful to obtain analytical results on the 
Efimovian trimer spectrum in the large / limit, see sub- 
sections HV^ and [TVBl 



B. Zero range case (i?* = 0) at negative energy 

Once the imaginary values of the Efimov exponent s 
are determined by solution of the transcendental equa- 
tion A/(s) = 0, which is possible for I odd and a larger 

than a critical value ai , see subsection llV Al one can de- 
termine the corresponding Efimovian trimer solutions of 
the zero-range theory at arbitrary energies, in free space 
and also in an isotropic harmonic trap, using the general 
real space formalism relying on separability of the Bethe- 
Peierls problem in hyperspherical coordinates |27l. [29l . [30| . 
Here we find it interesting to deduce from the real space 
solution the explicit form of the momentum space solu- 
tion, restricting for simplicity to the values of the mass 
ratio a and the (necessarily odd) angular momentum / 
such that the Efimov effect takes place. This may be 
useful for example to calculate the atomic momentum 
distribution of the Efimov trimer states in the zero-range 
limit, as was done for three bosons in [3J]. 

The idea to obtain DCk) is to take the limit of the 
position ri of the extra particle and the position r2 of 
a fermionic atom converging to the same location, with 
a fixed value of their center of mass position and of 
the position r^ of the second fermionic atom. Accord- 
ing to the Bethe-Peierls framework, the atomic wave- 
function ■(/;(ri,r2,r3) shall then diverge as 1/r, with 
r = ri — r2, with a factor depending on the Jacobi coor- 
dinate X = r3 — (ilfri + mY2)/{M + m): 

X fixed '^[p^) 



V'(ri,r2,r3) 



(47) 



Then one calculates -4(x) in two different ways. First, 
one uses Efimov solution exposed in the Appendix O 
Prom Efimov's ansatz (jAll) and from the solution F{R) oc 
KsKm/iiY^^qR] of the hyperradial equation (|A6[) . where 
rh is an arbitrary mass unit and Ks a Besscl function, 
one finds for a vanishing angular momentum along z: 



where we have set 



1/2 



(l + 2a)i/2 



q = 



(48) 



(49) 



where /Zam is the reduced mass of one fermionic atom 
of mass m and one pair of fermion-plus-extra-particle of 
mass m + M. Second, after inspection of (fT3| . one cal- 
culates the atomic wavefunction -0 by taking the Pourier 
transform of A(ki, k2, k3)(27r)^5(ki -I- k2 + k3), where A 
is given by of (|16|) , taking the functions x equal to unity. 
As expected, the contribution involving x(ki2)I?(k3) is 
the one leading to a divergence of ip for r — > 0. Inte- 
gration over ki is straightforward thanks to the factor 
(27r)'^(5(ki -I- k2 -I- ks). Then integration over k2 can be 
done after the change of variable k2 = ka — j^ks, using 



/ 



d^fc' 



-iU-^-r 



'Kr 



(27r)3 k'^ + 10 A-KT 

where K = [q^ + ^^^klY'^ > 0. Por r ^ 
approximates e^^'"/r with 1/r and one obtains 



A{^)^j^,Di).,y-^-. 



(50) 



0, one 



(51) 



We thus set, for an arbitrary choice of normalization lead- 
ing to a dimensionless function: 

Z?(k) = /'d3xe-*--rAx)^-^^J^. (52) 

A first technique to calculate the integral in (|52|) is to 
use the expansion of the plane wave on spherical har- 
monics: According to the identity (8.7-13) in [23|, the 
amplitude of the function x — )> e^^'^ on the spherical har- 
monics Y;'"'(x) is 47ri'ji(fea;)[y;™'(k)]*, where the spheri- 
cal Bessel function is real and may be expressed in terms 
of the usual Bessel function J: 



^'^^'^^^-^V-^'^"'^^' 



(53) 



It turns out that the integral over M+ of the product of a 
power law and of two Bessel functions may be expressed 
exactly in terms of the hypergeometric function 2^1, see 
relation 6.576(3) in [3l| : 



i?(k) = 2^3/2(-z)'(fc/Q)'y;"(k) 



PI 



^+s^ 



1 + s 

y-2Fi i + — — ,1 



l-s , 3 

2 ' 2 



P(; + 3/2) 
fcVQ' 



(54) 



This immediately shows that D{k) vanishes as /c' for k -^ 
0, which is generically the case for a regular function of 
angular momentum I. It also allows to obtain the large 
momentum behavior of I?(k) [35i] , 



/?(k) 



2n2 



ilT^Q 



k/Q- 



(l 




fc2 cos(7rs/2) 

{l-l)/2 



Y,°{k) 



n 

n=0 



s- (2n+l) 



2n 



m 



4(l-5)(fc/Q)2 



oiQ/kr 



+ c.c^ (55) 



which will play a crucial role in what follows to obtain 
the three-body parameter at non-zero i?*. In particular, 
Eq. (f55|) shows that A:^Z3(k) is asymptotically an oscillat- 
ing function of k/Q that is log-periodic: The same pat- 
tern is reproduced when k is multiplied by exp(27r/|s|). 
For fc^_D(k) to approach this asymptotic oscillating func- 
tion, the very stringent condition k/Q ^ exp(27r/|s|) 
is fortunately not required, it is simply sufficient that 
k/Q ^ 1 (for I and \s\ not much larger than unity). 

A second technique to calculate (15^ is to use spherical 
coordinates of axis the quantization axis z. The integral 
over the azimuthal angle (f> is straightforward. The in- 
tegral over the modulus x can be performed if one uses 
the integral representation 8.432(1) of the Bessel function 
given in [31|, Ks{z) = Jq °° dt exp(— zcosht) cosh(s<). In 
the integral over the polar angle 6, one replaces i^;°(x) by 
its expression in terms of the Legendre polynomial Pi (u) 
with the variable u — cos 9, and one integrates by parts 
the factor l/(coshi + iku/Q)"^ that appeared after inte- 
gration over X. The integral over t can then be performed 
with contour integration: 



+ 00 



dt 



-iSt 



_ sin[S'(/3 -I- i7r/2)] 
2tt cosh t + i sinh /3 i sinh(S'7r) cosh /3 



(56) 



where S and /3 are real quantities. Changing to the vari- 
able (3 = a.vgs\i(ku/Q) in the integral over m, one obtains 



i?(k) = y,°(k)- 



2i7r2Q2 



cos[S'argsh(fc/(5)] 
fccosh(S'7r/2) I (fc2 + Q2)i/2 
1 ^aigsh(fe/Q) /r)sinh/?\ 1 



(57) 



where we have set s — iS. Generalizing the technique 
of subsection IIII Al to the derivative P{ of the Legendre 
polynomial allows a direct evaluation of £'(k) without 
the need of hypergeometric functions. 



C. Case _R» > at zero energy 

In this subsection, as in the previous one, we restrict to 
the case where an Efimov effect takes place: Anticipating 
on results of subsection IIV A[ the angular momentum 
quantum number I is odd and the mass ratio a is larger 
than the corresponding critical value ac\ so that the 
function A;(s) has a single purely imaginary root si — iSi 
with positive imaginary part: 

Ai{iSi) = with Si > 0. (58) 

As remarkably shown for three bosons in 12!] the narrow 
Feshbach resonance model can be solved analytically at 
zero energy, which gives access to the three-body param- 
eter, that is to the energy scale in the asymptotically 
geometric Efimovian trimer spectrum. 

The underlying idea of the solution is that, at zero 
energy and after division of the overall Eq. ([M|) by a 



factor k, the integral part of the resulting equation is 
strictly scaling invariant. If one takes as variable the 
logarithm of k rather than k, this strict scaling invariance 
corresponds to a translation invariance, which suggests 
that the integral part is simply a convolution product and 
leads one to perform a Fourier transform with respect to 
Infc. More precisely, we use the ansatz: 

/«(fc) =e-2='F(')(a;) with x = In (fci?, cos i/) , (59) 



where the factor 



shall lead to a convolution kernel 



with the desired even function, and cos i' is a function 
of the mass ratio given by Eq. PO)) . Dividing Eq. ([M)) 
for i? = by fc, injecting the ansatz ((59|) in the resulting 
equation, and finally multiplying by e^^, we obtain 

/ + 00 
dX ICiiX - x)F^'\x) 
-oo 

(60) 
with the kernel 



Kiix)^lfdu- ^'(") 



27r J_i cosh(x) + jt:^u 



(61) 



As expected, (|60|) is the convolution product. We thus 
introduce the Fourier representation of the function F^''' 
M: 



F^'\a 



+ oo+iO+ 



iO+ 



dS 

( 

27r 



,^Sxp{l)(^gy 



(62) 



We expect that the function F^'^(a;) oscillates periodi- 
cally for X -^ — cxd: When x — >■ — oo, the momentum 
k tends to 0, it becomes much smaller than 1/-R»; one 
enters a universal zero-energy, zero-range regime where 
solutions of the type (|35t are obtained, with s purely 
imaginary for an Efimovian solution. According to (|59p . 
this implies that _F(')(x) has plane wave oscillations at 
X -^ — cxD, with a wavenumber ±Si since both s = iSi 
and s = —iSi are roots of A/(s) = 0. More precisely, we 
expect that there exist coefficients A± such that 



F«(a 



^iSix 



A, 



-iSiX 



0(1). 



(63) 



This will be checked a posteriori. As a consequence, the 
Fourier transform F^^'>{S) has singularities on the real 
axis. More precisely, Eq. (l63l) leads to the natural con- 
clusion that 



F'^^^{S) has simple poles in S ^ ±Si. 



(64) 



This is why the integration contour in (j62p is infinites- 
imally shifted upwards in the complex plane. This 
is a standard procedure in physics, see for example 
the expression of the unitary evolution operator as a 
Fourier transform of the resolvent for a system with time- 
independent Hamiltonian [37| . In this Fourier representa- 
tion, the convolution becomes a product, and one needs 
to calculate the Fourier transform ICi{S) of the kernel 



function K. Exchanging the integration over u and x, 
one then recovers exactly the integral in Eq. ([57]) know- 
ing that s — iS in that equation, so that 



, Im z 



K:i{S) =A{iS)-coav. 



(65) 



The only subtle part is the determination of the Fourier 
representation of the function x -^ e^F^^'>{x). Multiply- 
ing (|62|) by e^ gives 



2"i^(')(a;) 



+oo+iO"f 



D+iO+ 



dS 
2^' 



i{S-i)xp{l)/g\ 



(66) 



This is not directly of the Fourier form (|5^ because S — i 
rather than S appears inside the exponential. Now, if 
the integrand, that is here in practice the function S -^ 
F^'''{S), is a meromorphic function with no singularities 
(no poles) in the band 0+ < Imz < 1 + 0+, that is 
< Imz < 1, where z £ C, we can shift the integration 
contour in (|66p upwards by one unity along the vertical 
axis in the complex plane, so as to replace S — i with S. 
Under the hypothesis 

F^ -'(z) has no singularities for < Imz < 1, (67) 



we thus have 



X 



+oo+iO+ 



oo+iO+ 



dS 

2^ 



e^Sxpii)^g^-^^ 



(68) 



which is exactly of the Fourier functional form (|62|) . that 
is the function S —¥ F^'^ {S + i) is the Fourier transform 
of the function x — > e^F^^>{x). This procedure is sum- 
marized on Fig[T] where the original ((66|) and the shifted 
(j68| integration contours are plotted, and where the lo- 
cations of the poles of -F'''"(z) are also indicated. Eq. (|60|) 
then reduces to 



= f(') {S + i) cos v + Ki {iS)F'^^^ (S) 



(69) 



for all real S. 

To solve §^ we follow [12.]. 
ansatz 

F<^^\S) = 



We first introduce the 



sinh[7r(S' + 5*0] 



Ci{S) 



(70) 



where we recall that iSi is the positive-imaginary-part 
root of Ai{s). The factor with the hyperbolic sine is 
carefully chosen so as to give a minus sinus under trans- 
lation S ^- S + i, a priori introducing poles in z = —Si, 
z = —Si ± i, z = —Si ± 2i, etc, for the function 
F*^'^ (z). This does not introduce singularities in the band 
< Imz < 1 provided that the function Ci{z) tends to 
zero for z ^ —Si + i, a point to be checked a posteriori. 
The unknown function Ci solves 



Ci{S + i) cos iy = Ai{iS)Ci{S) 



(71) 



on the real axis. According to the expectation (|M)) . the 
function Ci{S) has no pole in 5' = —Si, since the 1/sinh 



original contour z=S+i 



-S, 



-1/2 



shifted contour z=S+i+i 

i I 



1/2 



Rez 



FIG. 1: (color online) In the complex plane, illustration of 
the procedure used to determine the Fourier transform of the 
function x — > e^F"'{x): Under the a posteriori checkable 
hypothesis (|67p (no pole of F^ \z) in the grey area), one can 
shift the integration contour in (I66|) upwards by one unity 
along the vertical axis, to obtain the Fourier representation 
Eq. dMl) ofs -^ e="F(')(s). Circles: Poles of F*'' (z) due to the 
function hyperbolic sine in the denominator of (1701) . Pluses: 
Poles of F^''\z) due to the factors r(l + iz — iv„), n > 0, 
in the numerators of the infinite product (|73p . Black crosses: 
Poles of F^'\z) due to the factor r{iSi - iz) in ((73|. Colored 
crosses on the Imz < part of the imaginary axis: Poles of 
F'^'\z) due to the factors r(— m,i — iz) for n > in (f75)) [red: 
n — 1; blue: n — 2; violet: n = 3; cyan: n = 4]; from top 
to bottom, n — 1 (three crosses), n = 2 (two crosses), two 
triplets n = 1, 2, 3 and a quadruplet n = 1, 2, 3, 4. The figure 
corresponds to I = 1 and a = 20. 



factor in ((70| already has a simple pole in S = —Si, 
and Ci{S) has a simple pole in S = Si, since the 1/sinh 
factor has no pole there i38i] • As in |12j] one then uses the 
Weierstrass representation 



H^s) 



COSJ^ 



n 

neN 



S^-vl 



(72) 



where the overall factor cos;^ is the limit of A/(iS') for 
S ^ oo, see Eq. (gH). In Eq. ([71]), Uq ^ Si is the real 
positive root of the function 5* — ;• Ki(iS), —uq is the 
real negative root, and the ±u„, n > 1, are the purely 
imaginary roots of 5" — >■ Ai{iS), ImM„ > 0. The w„'s 
are sorted by ascending order of their imaginary parts, 
and for general values of the mass ratio a, they form an 
irregular, aperiodic sequence. They of course depend on 
the angular momentum /. On the contrary, for general 
values of the mass ratio a, the poles of S* ^' A{iS) are 
found from (|151) to sim ply be ±w„, with w„ = i{2n + l + l) 
for all integers n > [39|. Finally, one can check as in 
[13 | that the function Ci is given by the infinite product 



C'liS) = 



r{iSi-iS)T{l + iS-ivo) 
r{l + iS + iSi)T{-iS - ivo) 

n T{-iS - iUn)T{l + iS - ivn) , , 



neN 



^ T{—iS — iVn)T{l + iS — iu^) ' 



where we recall that the m„ and Vn depend on L In 
particular, one can check that this expression vanishes 
for S* — > — S*; + i, as required above Eq. (|7T|) . and has no 
pole in the band < Im z < 1 of the complex plane [40] ■ 
Together with ((70|) . this constitutes the desired solution 
at zero energy for i?* > 0. 

An important application of this result is to calculate 
the previously mentioned low-fc or cc — >■ — oo behavior 
(|55| of the solution [42l , which is a universal regime that 
has to match the zero-range model. Since x < 0, in 
applying the usual contour integration technique to ()62p , 
we close the integration contour following a half-circle in 
the lower part Im z < of the complex plane. According 
to the Cauchy residue formula, one gets for F'^^^x) a sum 
of terms proportional to e*^"^, where the sum is taken 
over all poles z„ of the integrand in the lower half plane 
(see the pole locations in FigH]). For a; — > — oo, the poles 
with a non-zero imaginary part have a contribution that 
vanishes as O(e^), and the sum is dominated by the two 
poles z — ±Si on the real axis, which are the only ones to 
give purely oscillating, non-decaying contributions. The 
corresponding residues of F'^^^S) can be deduced from 



F«(5) 



Cii-Si) 
sS'-Si S + Si 



and F^-^^{S) 



[Ci{~Si)Y 

s^s, S - Si 



(74) 
The value of the first residue directly results from the 
ansatz ((70)) and the absence of pole of the function Ci{S) 
in S = -Si. If one further uses dnl) for S ^ Si, 
one finds for the second residue [7r/sinh(27r5;)]Ci(S'; -I- 
i) cos v/ih'i{iSi). Properties of the explicit form (1751) and 
of the Gamma function, as in [12| , lead to ([74| |43| . Fi- 
nally, turning back to the k variable and to the function 
D: 



D(\i) = ~i (-^ 



Yn^) 



X {[Cii~Si)]*iqR,k/Qy'^'+c.c. + 0{qR,k/Q)} (75) 



where we recall that Si > and q/Q = cosi' as in (j49| . 
The asymptotic form in the right-hand side of ([75]) is sat- 
isfactory: It is indeed a superposition of solutions of the 
zero-range model at zero energy, see (|35|) . A first impor- 
tant point is that it is actually a specific linear combina- 
tion of the solutions with exponents s — iiS*;, with rel- 
ative amplitudes depending on the Feshbach length _R*. 
This selection of the right linear combination amounts to 
adjusting the three-body parameter to its right value in 
the Danilov three-body contact conditions 28] . Eq. ([75]) 
thus constitutes a microscopic derivation of this three- 
body parameter in the limit of an infinitely narrow Fes- 
hbach resonance fill, lOJ. A second important point 
is that the solution -D(k) starts approaching the log- 
periodic oscillatory asymptotic form (fci?* cos:^)^*'^' as 
soon as kR^ cos v < 1, there is no need to require that 
fci?, cos ly < e^^^^^' : There is no need to require that the 
oscillatory form has performed at least one oscillation to 
have I?(k) well approximated by it. 



IV. ANALYTICAL RESULTS ON TRIMER 

STATES 

In this central section of the paper, we develop a phys- 
ical application of the particular analytical solutions of 
the previous section. From the zero-energy and zero- 
range solution, we first determine, for each value of the 
angular momentum I, the critical mass ratio ac lead- 
ing to the Efiniov effect and the corresponding purely 
imaginary Efimov exponent si = iSi, with the conven- 
tion Si > 0. Accurate asymptotic estimates for these 
quantities are obtained and are compared to the Born- 
Oppenheimer approximation. Second, from a matching 
of the {E < 0, i?, = 0) solution to the (^ = 0, i?* > 0) so- 
lution, we determine the Efimovian part of the spectrum, 
in particular the global energy scale ^iiobai appearing in 
that geometric spectrum. The dependence of -Egio^ai °^ 

the mass ratio, close to the Efimov threshold a -^ ac 
and arbitrarily far from it {a -^ +oo), is studied. Also 
the variation of -Egio^ai with the angular momentum I, in 

particular for a close to the critical mass ratio ac , is an- 
alyzed. Third, using the Born-Oppenheimer approxima- 
tion expected to be asymptotically exact for a diverging 
mass ratio, we show that the hydrogenoid character grad- 
ually takes over the Efimovian character in that limit, ex- 
cept in a vicinity of the E = accumulation point (which 
remains Efimovian). 



A. Efimovian threshold ai and exponent si 



As a physical application of Eqs. (|43I44I46I) . we deter- 
mine the values of the mass ratio a and of the angular 
momentum I such that the Efimov effect takes place, that 
is the transcendental equation A;(s) — admits some 
purely imaginary solutions. 

A useful guide is the Born-Oppenheimer approxima- 
tion [5| that becomes exact in the limiting cases of van- 
ishing or diverging mass ratio a. It indicates that the 
Efimov effect should take place for large enough values 
of a and for odd values of L In the limit a — J> 0, the 
extra particle is indeed infinitely massive, so that the 
fermions see a fixed point-like scatterer with infinite scat- 
tering length, which does not support bound states. In 
the opposite limit a — )■ +oo, the extra particle sees the 
very massive fermions as two fixed point-like scattcrcrs 
of positions r2 and r^ , with which it forms a single bound 
state, of energy 



eo(?'23) 



'2M^ 



(76) 



in the zero-range Bethe-Peierls model, and with a wave- 
function that is symmetric under the exchange of r2 and 
r^. Here the constant C obeys 

C = exp(-C) so that C = 0.567 143 290 409 .. . (77) 
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It can be related to the Lambert function by C = VF(1), 
and it is sometimes called the J7 constant. Since the 
global state vector is fermionic, this extra particle wave- 
function can be combined in the Born-Oppenheimer fac- 
torized form with an odd orbital fermionic part only. The 
effective potential seen by the fermions is thus the sum 
of the Born-Oppenheimer potential eo(r23) and of the 
angular momentum centrifugal part h?l(l + l)/('7ir23). 
In three dimensions, the zero-energy solution in that ef- 



fective potential has to be written as 



23 



-1/2 



to match 



the usual definition of the Efimov exponent which is 
given in two dimensions [44| . We then obtain the Born- 
Oppenheimer approximation for the root s: 



„2 

sbo 



1 + 



1 



1 



aC^ {I odd) 



(78) 



to be used in the regime where s^q < 0. Interestingly, 
this Born-Oppenheimer approximation can also be used 
for a non-zero R^ to obtain exact results on the non- 
Efimovian low-energy trimers for a -^ -l-oo, as developed 
in subsection IIV CI 

Turning back to the exact equation A/(s) = 0: For a 
fixed value of I, we expect that the solutions s of the 
transcendental equation are continuous functions of the 
mass ratio a. The critical values ac of a for the emer- 
gence of the Efimov effect are thus such that A;(0) = 0. 
From a numerical calculation of Ai{s = 0) as a function 
of a ranging from to -l-oo, we indeed find, for even Z, 
that A;(s = 0) has a constant positive sign which means 
the absence of Efimov effect. For each odd /, we find 
that Ai{s = 0) changes sign once (from positive to nega- 
tive for increasing a). The resulting values of ai are for 
example: 



a^'-^) =13.60696..., a 
a('=5' = 187.9583 . . . , a 



(1=3) 
(i=7) 



75.99449 . . 
: 349.6384 . . 



(79) 



The critical mass ratio for I = 1 was already given in 
Q, and in 45| for I = 3 and I = 5. For larger / we 

have checked that [ac C'-^/2]^/^ is indeed very close to 
the approximation I + 1/2 resulting from ([75)1 . For these 
few odd values oil, we then solve numerically A; (iS) — 0, 
where 5' > 0, to obtain the Efimov exponent as a function 
of a for a > ac . For each values of I and a > ac , it is 
observed that Ai{iS) for S* > is an increasing function 
of S, so that a single pair of iiS*; imaginary roots is 
obtained. The results are shown as solid lines in Figl2] 
For comparison, the Born-Oppenheimer approximation 
([75]) is plotted as dashed lines in that figure. As expected 
intuitively, it approaches the exact result in the large I 
(and thus large a) limit. 

One can however be more precise in the evaluation of 
the accuracy of Eq. ((75)) . For a fixed value si = iSi of the 
Efimov exponent, one can perform a large I expansion of 
the mass ratio, as shown in Appendix iDl to obtain: 



\c-'ia 



v(0^ 



Sf [I + 0(1/1)] 



(80) 



where the critical mass ratio ai has the large I expansion 
lc'ai^\^^(^l + ^y + A + 0{1/1). (81) 

A is thus the first correction to the critical mass ratio 
predicted by the Born-Oppenheimer approximation (|78p . 
Accidentally it has a very small numerical value |46|: 



17 -g^ 
12 



^- 



1 



C + 1 



-0.016 259 165. 



(82) 

so that the Born-Oppenheimer approximation for ac is 
in practice quite good for large I. For a — > +oo for a 
fixed angular momentum Z, one can also obtain from the 
results of Appendix[E]the following asymptotic expansion 
for the imaginary root iSi of A/(s): 



S? = 



C^ 



„ a- n+ - 

Q— >--i-oo 2 \ 2 



^ + 0(l/a) 



(83) 



where the first two terms in the right-hand side constitute 
the Born-Oppenheimer approximation and the angular- 
momentum independent constant 5 is different from zero: 



12 



4 C^ 



31 + C 



- = 0.657 684. 
4 



(84) 



Note that Eq. (1551) is not in contradiction with Eq. ([50)) 
[combined with Eq. ()81l) ] because the considered limits 
are totally different, Si/l either diverges or vanishes. As 
s^ is a coefficient in the Born-Oppenheimer potential, 
a suggestive picture obtained from Eq. (1551) is that the 
Born-Oppenheimer approximation has an error 0(oP) on 
the effective potential seen by the fermions. Mathemat- 
ically, Eq. (|55)) shows that the error on Si due to the 
Born-Oppenheimer approximation vanishes as 1/a^", 
and that the approximation gets the correct leading term 
a^' ^ and the correct subleading term a" (which turns out 
to vanish) in an asymptotic expansion of Si in powers of 
a^/^. These remarks are useful for subsection IIV CI 



B. Efimovian part of the trimer spectrum 

In presence of the Efimov effect, we have obtained so 
far two solutions to the integral equation p4p in limit- 
ing cases, see subsection IIIIBI for i?* = 0, i? < 0, and 
subsection llll CI for i?* > 0, i? = 0. How can we then ob- 
tain an approximation for the corresponding spectrum of 
trimers, which requires to have both i?* > and E < 07 

Roughly speaking, dropping for simplicity a possible 
dependence on the mass ratio a of the various bounds, 
the solution of subsection IIIIBI is expected to constitute 
an accurate approximation of the trimer solution at in- 
terparticle distances much larger than R^ , that is at mo- 
menta smaller than h/R^. In a symmetric manner, the 
solution of subsection IIIICI is expected to well approx- 
imate the trimer solution at short enough interparticle 
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range of applicability of 
solution E=-(!iq) /2|a., R^=0 







its asymptotic range : 
:< > 



matching interval 



1/R. 



its asymptotic range 



range of applicability of 
solution E=0, R,,>0 



k 

CX) 



FIG. 2: Values of the Efimov exponent s as a function of the 
fermion-to-extra-particle mass ratio a — m/M, for several 
values of the angular momentum I. The modulus of s is shown 
only over the interval of mass ratio where s is purely imaginary 
(that is where the Efimov effect takes place). This was found 
to occur only for odd values of I, and for a single pair of ±s 
roots of Ai(s) = 0. Solid line: Numerical result from the 
exact function A;(s) as given by (|44p . Dashed line: Born- 
Oppenheimer approximation (|78|l . 



distances, smaller than 1/q, where the E < wavefunc- 
tion only weakly departs from the zero-energy one. This 
corresponds to momenta much larger than hq. There thus 
exists an interval of momentum hk over which both lim- 
iting solutions are close to the physical trimer solution, 
q < fc < 1/R*, if 



qR* < 1. 



(85) 



In this case, we can match the two limiting solutions, 
as summarized in Fig. [31 Over this matching interval 
of momentum, one has A; ;^ g so that the solution of 
subsection nil Bl is in its large k regime given by (|55p with 
s = iSi. Over this interval, one also has k ^ 1/^*, so 
that the solution of subsection IIII Cj is in its low-fc regime 
given by (j75p . These two limiting regimes are compatible 
(within an arbitrary normalization factor) if q is of the 
form El 



In — yglobal^ 



''/ki 



n integer. 



(86) 



where we recall that the Efimov exponent is si — iSi, 
Si > and the trimer energy is E = —h'^q'^/{2fi), with 
/i = mM/{m + M). Eq. ((86)) is the expected geometric 
Efimovian spectrum, and we have the explicit expres- 
sions for the corresponding global wavenumber and en- 
ergy scales: 



9, 



(0 



■p^'/l"'! and K^'^ - 



2^2 



^global ^^ - — --global ^^2 

where 9i is the phase of the complex number 



o2e,/|s,| 



(87) 



Zi 



\Zi\e 



i0l 



n("-^') 



siCii-Si) (88) 



FIG. 3: (color online) Matching of two limiting solutions of 
the integral equation (I34|l . the solution {E = —h q /{2fi) < 
0, -R* = 0) [see text in black, in the upper half of the figure] 
given by Eq. (f51)) and the solution (E = 0, R, > 0) [see 
text in red, in the lower half of the figure] deducible from 
Eqs. (|29I59I62I70I73|I . over a common interval of values of k 
[the "matching interval" , blue segment on k axis] where they 
are both in their asymptotic regimes Eq. (|55p and Eq. (|75|) . 
This matching procedure leads to the Efimovian spectrum 
formula (|86p with a global scale given by (|87p . This procedure 
makes sense when g_R* <^ 1, and it is expected to be exact 
when qR, — >■ 0, that is for the quantum number n tending to 
infinity for a fixed (purely imaginary) Efimov exponent s; (as 
in Efimov's historical solution) or (less usually) for |s;[ tending 
to zero at fixed quantum number n > 1. For simplicity of 
the figure, we have dropped factors slowly depending on the 
mass ratio a, such as cos v, that is we have assumed that the 
angular momentum I and |s;| are not much larger than unity. 



and the function Ci{S) is given by Eq. ([73]) . At fixed 
angular momentum, the phase 9i depends on the mass 
ratio a — m/M . For a tending from above to the critical 

value ac , Si —>■ 0, we show in the Appendix IBJ that Zi 
tends to a real and positive number. We thus choose the 
usual determination 9i — Arg Zi in that limit, and extend 
it by continuity to all larger values of a. 

For the analytical developments that follow, obtained 
for the global scale ggio^ai i^ ^^^ limit of large angular 

momenta for a — ac fixed, or in the limit of a large mass 
ratio at fixed angular momentum, there is an actually 
more operational expression for 9i , that does not require 
the determination of all the roots and poles of Ai{iS) to 
evaluate Ci{~Si). As shown in Appendix [Cl one has the 
series representation of the phase 9i of Zi : 

9i = Im[ln r(l + ^5*0 +hi r(l + 2iSi) + 2 In r(/ + 1 - iSi) 



fc>i 



lnr(/ + 2-iS'/)] + 



Si 

dS In 



Ai{iS) S^ + il + iy 



(2fc)! d52'=-i 



In 



cosiy S^-Sf 
AiitS) S^ + {l + lf 



cos^ 



S^-Sf 



S=Si 

(89) 
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where we recall that Si is the positive root of the func- 
tion Ai{iS) and the B^ are Bernoulli's numbers, Bi — 
— 1/2, B2 — 1/6, . . .. Since (|89|) was obtained from Stir- 
ling's series, which is an asymptotic series, we expect that 
((89)) is also an asymptotic series. We thus investigated 
numerically how many terms one has to keep in practice 
to have a good accuracy. In FigH^, it is found remark- 
ably that the zeroth-order approximation, consisting in 
omitting all the terms in the sum over k in (j89p , already 
gives in practice a sufficiently accurate approximation for 

the global scale qgiohni- ^^ FigHb, it is shown that the 
difference between the exact value of 9i, obtained from 
the infinite product representation of Ci{—Si), and the 
zeroth-order approximation, omitting all k terms in (1891) , 
is non-zero but is very accurately accounted for by the 
term fc = 1 in ((89)) . 

From the zero-range theory by Efimov, it is expected 
that the geometric form Eq. ([55)) of the spectrum is 
asymptotically exact in the limit of a large quantum num- 
ber, n — s> +00, see the constraint (|85)) . How large the 
values of n should be to reach this geometric behavior 
will be evaluated numerically in section |V] [in the mean 
time see the discussion that follows Eq. (IHT)) ] . The advan- 
tage of the non-zero i?* calculation is that it also gives 
the global energy scale in the spectrum, see the factor 

^global i'^ (|55)) explicitly given by (|57)). As compared to 
the bosonic case [il|,[l4|: there is an additional knob here, 
which is the mass ratio a. How does the global energy 
scale vary with the mass ratio ? 

For a — ?• ac ■ The behavior of the global energy scale 
close to the critical mass ratio can be determined from the 
results of Appendix [B] where it is shown that 0i vanishes 
linearly with Si: 

n 

lim — ^ = -ij{l + 1) + iil){l) - i;i~ivo) - V'(l - ivo) 

a^a<" \si\ 

+ 00 

n=l 

(90) 

Here u'^,'n > 0, is the value of the complex root m„ of 
the function S — > Ai{iS) at the critical mass ratio, the 
Vn = i{2n -|- ^ -|- 1), n > 0, are the poles of S" — > Ai{iS), 
and the function il^^z) is the digamma function, that is 
the logarithmic derivative of the F function. As a conse- 
quence, the global energy scale has a finite limit at the 
threshold for the Efimov effect! We give here a few cor- 
responding values obtained from the rapidly converging 
formula (IB6I): 



,0=1) 



1-2 „('=3) 



9giobai^* ^ 6.56577 • 10-^ q^-^R* ^ 6.12349 • 10"^ 
&iR* ^ 1-62809 • 10-3, g(;=J)^i?, c 6.48952 • 10"^ 



(91) 



At this threshold, |si| — ;> so we expect that Eq, 
becomes actually exact for the quantum number n > 1 



since gn i?* tends to zero in that limit, whereas (|86)) is 
clearly invalid for n < —1. The case n = is dubious: 
Although Qq R^ does not tend to zero at the threshold, 
it assumes very small values, so may be Eq. ()86l) still 
makes sense. We can not however say more at this stage, 
and the question whether the quantum number n ~ 1 
corresponds or not to the ground trimer state (for a given 
angular momentum I) will be answered in section IV) 
For / — > -|-oo: Another interesting question is to deter- 
mine how (/global depends on the angular momentum I 
at a fixed distance of the mass ratio a from the critical 
value ac , that is roughly at constant values of the Efi- 
mov exponent si. After a numerical evaluation of ()87p . 

as detailed in the note [33], we found that ggio^ai drops 

rapidly for increasing I, roughly as l/(ac )'^^^- Accord- 
ing to ([75)) the critical mass ratio scales approximately as 
(l + 1/2)^, so that we expect that q^^^g^i approximately 
scales as l/{l + 1/2)'^, an approximation that becomes 
rapidly excellent with increasing I as soon as I exceeds 
unity, see FigO 

The scaling of gl-iobai as l/l"^ at a fixed distance from 
the critical mass ratio can be obtained analytically for 
I — > +CXD using the expression ([5^ for the angle 9i and 
the expression Eq. (I46p for Ai{iS), as detailed in the Ap- 
pendix [d) 



^global-"-* 



j-a« fixed (1 + C) e^^ 



,,,^^^-^i^-ni + ^SiHinrii + 2^Si)y^^ (92) 



where 7 = 0.577 215 664 9 . . . is Euler's constant and 
si — iSi, Si > 0. In that limit, one can use the Born- 

Oppenheimer-type relation, Sf = {a — ac )C'^/'2,, as 
shown by Eq. ()80p . The asymptotic result (1^^ is plotted 
as a dashed line in Fig(5] and well reproduces the large-Z 
numerical results. 

For a -^ -t-00: The behavior of the global energy scale 
in the limit of an infinite mass ratio (for a given I) is 

determined in the Appendix [Ej It is found that q„ia^,gx 
has a finite limit, which remarkably is also independent 
of the angular momentum I : 



,« 



global-"-* 



i?* 



a— >-+oo 



2(1 + C)e" 



3.31582. 



(93) 



where C is defined by ([77)) and J is the integral 



J 



c 



dx 



1 1 



1 



0.05630577. 



C l-xe^ C -X 

(94) 
As we shall see in subsection II V CI for a — >■ -|-cx3, the low- 
lying part of the trimer spectrum is hydrogenoid rather 
than Efimovian, so that (1931) is relevant only for trimers 
with diverging quantum number n. 
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FIG. 4: (color online) In (a) value of the global wavenumber scale qJ^-^^T^ of the Efimovian part of the spectrum, see (|87|l . as 
a function of the difference of the mass ratio a from its critical value, for all odd angular momenta up to Z = 7. Thin solid 
lines: Exact result obtained from (|88[) and from the infinite product representation (|73|l [the more rapidly converging form 
(|C3IC5|) was used in the numerics]. Thick dashed lines: Zeroth-order approximation omitting all the terms in the sum over 
k in (|89p. In (b) the (small) difference between the exact 6i and the zeroth-order approximation Q\^^ is plotted as thin solid 
lines, and the value Qi "^ of the fc = 1 term of (|89p is plotted as thick dashed lines. In (c) the values of giiobai resulting from the 

Born-Oppenheimer plus semi-classical approximation are given in dashed lines [see Eq. (|108p ]. as functions of l/(a — Oc )^", 
showing how the /-independent q — > +oo limit is reached. The solid lines correspond to the exact numerical data of (a) and 
the dotted straight lines corresponds to the asymptotic expansion Eq. piO[) . 




a-a 



FIG. 5: Dependence, with the angular momentum /, of the 
global wavenumber scale qJ^-^^T^ in the Efimov trimer spec- 
trum, as obtained from the exact result (|87p . The figure shows 
the quantity (f^ohaJ^* multiplied by (/ -I- 1/2)^, as a function 

of a — Oc , for (odd) angular momenta equal to / = 1, Z = 3, 
/ = 5 and 1 = 7 from bottom to top. This quantity is ob- 
served to converge rapidly to a finite limit for increasing I. 
The dashed line represents the analytical prediction (I92p for 
that limit. We recall that £giobai = —^Q^oha.i/i'^iA- 



interaction potential on a narrow Feshbach resonance. 
Then from Schrodinger's equation for the relative 
wavefunction of the two fermions, 



E^{V2Z) = Ar,3V'(r23) + e(r23)V'(r23) (95) 

m 

which is in the odd I sector due to fermionic antisym- 
metry, we determine the low-energy trimer states in the 
limit a -^ -|-cx). 

We first fix the positions of the fermions to V2 and r^ . 
Each fermion then acts on the extra particle as a fixed 
scatterer of infinite scattering length, zero true range but 
finite effective range r^ = —2R^. The effect of such a 
scatterer is then represented by modified contact condi- 
tions on the wavefunction (j){ri) of the extra particle in 
the so-called effective range approach 7, 48] . For a bound 
state 0(ri) of eigenenergy 



e(r23) 



2^2 



' 2M 



(96) 



where the dependence of k > on r23 is for simplicity 
omitted in the writing, we impose the boundary condi- 
tions 



Born-Oppenheimer approximation and 
hydrogenoid trimer spectrum 



As explained in subsection IIVAI the Born- 
Oppenheimer approach is a natural tool when the 
fermions become arbitrarily massive. The extra particle 
then mediates an attractive interaction potential e(r23) 
between the heavy fermions, with r23 = r2 — r3 the 
relative coordinates of the fermions. We calculate this 



(ri) = A,, 

ri —^Tn 



1 



n 



1 

flcff 



0(|ri-r„|) (97) 



in the vicinity of each scatterer n = 2,3. The effective 
scattering length is energy dependent, as l/oeS = 1/a — 
k'^re/2 for an incoming free wave of wavenumber k. Here 
the true scattering length a is infinite, the effective range 
is Tg = — 2i?* for a narrow Feshbach resonance, and k — 
IK is purely imaginary for a bound state, so that l/a^s = 
—K?Rf. In presence of the contact conditions, the extra- 
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particle wavefunction obeys Schrodinger's equation 



e(r23)0(ri) = - — 



3 

Ari(/)(ri) + ^47rA„5(ri 

n=2 



r„) 



mass me = 'm/2 and the CGS electron charge e such that 
e^ = h^/{2MR^). In terms of the wavenumber g intro- 



duced in Eq. (|26p we thus obtain the exact asymptotic 
result M 



(98) 
where the Dirac terms are due to the l/|ri — r„| diver- 
gences. The general solution is expressed in terms of the 
Green's function of the Laplacian at negative energy: 



Q: 



(Op 



a— ^+oo 



f- 

V8 



1/2 1 



n + l 



(104) 



(ri) = E^" 



-K|ri-r„| 



n=2 



|ri 



(99) 



Here the integer quantum number n start from n = 1 in 
each odd angular momentum sector Z, so that the hydro- 
gen spectrum reads — mee^/[2fi^(n + l)"^] with this con- 
vention. 

Since the eigenfunctions of the hydrogen atom are well 
known, it is possible to calculate the first correction to 

1/2 



The contact conditions ([57)) then impose (l + Ki?»)A2^3 = 

^3,2e"''''='V('«''23)- This 2x2 system has a non-zero so- the hydrogenoid spectrum, treating the l/r^^^ term in 

lution only for the symmetric case A2 = A^, where k 

solves 



1 + kR^ 



nr23 



(100) 



The wavefunction 0(ri) is then symmetric under the ex- 
change of r2 and r3. As the Born-Oppenheimer ansatz 
for the total wavefunction is ?/'(r23)'/'(ri; r2, rs), where 
the parametric dependence of (j) with the fermions po- 
sitions is made explicit, fermionic exchange symmetry 
indeed imposes V'(— r23) = — ■(/'(r23). 

For a fixed r23, Eq. (llOOp looks difficult to solve. How- 
ever, one can see that, for each positive k, it is solved by 
a single positive r23 [49|. Furthermore, rewriting k in the 
left hand side as (Kr23)/7'23, r23/R* may be expressed as 
an explicit function of u ~ Kr23, where u ranges from 
to the numerical constant C defined in (l77l): 



^^23 



(101) 



This allows a straightforward plot and study of the Born- 
Oppenheimer potential e(r23). We reach in particular the 
useful limiting cases 



e(''23; 



I'23->- + 00 






'^'^'^ r=^. -2Mm 



r23 



2R 



1/2 



1/2 

'^23 



0(1) 



(102) 



.(103) 



The asymptotic behavior (|102p reproduces the Born- 
Oppenheimer potential ([76| obtained for the usual 
Bethe-Peierls case R^ = 0. For a large enough a, and for 
each odd values of the angular momentum I, it ensures 
that the spectrum of Eq. (j95p is indeed Efimovian in the 
limit of large quantum number n, that is for _£ — > 0^. 
On the contrary, for fixed quantum numbers n and I, 
the bound states of ([95|l for increasing a = m/M are in- 
creasingly localized in the low r23 part Eq. (|103p of the 
Born-Oppenheimer potential. This means that the low- 
energy part of the spectrum is hydrogenoid, it becomes 
asymptotically equivalent for large a to the known spec- 
trum of the hydrogen atom if one takes for the electron 



Eq. (|103|) to first order in perturbation theory. The cal- 
culations are given in the Appendix |F1 we present here 
only the result: 



Q: 



" '"*- ~+oo n + l 



„ J Li 



1 ^(i)( ^i^ + l) \ 



1/2 



oil 

a 



?(0 



(105) 



where Rn is a rational number given by 






{n + 2iy. 



E 

fe=0 



(2fc)! 



2 A- 



(2A:-^1)2 



(fc!)2 

[2{2l + 1 + n - k)]\ 

{n- 1 - ky.{2i + 1 + n - ky: 



(106) 



An interesting question is to know whether Eq. (|105l) . 
that originates from the Born-Oppenheimer approxima- 
tion, is still exact |5l[. The discussion below Eq. ([M)) 
allows to hope so. We shall present numerical evidence 
in subsection IV CI that this is indeed the case. 

Another interesting aspect is to determine if the Born- 
Oppenheimer approximation, combined with a suitable 
semi-classical (WKB type) approximation, is able to 

determine exactly the global scale Qgiotai in the large 
mass ratio limit. For a fixed angular momentum /, 
using the technique presented in [52| and setting here 
E — —h^Q^/m, we obtain the semi-classical quantiza- 
tion condition [53| 



dr 



" 2/ \ (' 
2 ^ ' 



1/2)2 



-Q^ 



1/2 



= (n-l/2)7r 



where r„ii„ and r„iax 



(107) 
are the lower and upper roots of 
the integrand, and the quantum number n is any integer 
> 1. In the large n limit, Q tends to zero exponentially 
fast (this is the Efimovian part of the spectrum) , so that 
we can split the integral in (jl07p in two intervals, the 
interval rmin < ?" < ?'int where one can make the approx- 
imation Q ~ 0, and the interval rint < r < rmax where 
the Born-Oppenheimer potential can be replaced by its 
asymptotic expression (|102l) . The intermediate value rjnt 
simply has to satisfy R^ <C rjnt ^ |sbo|/Q where sbo is 
given by Eq. (|78p . Then the result does not depend on 
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the specific value of ri„t and one finds an approximation 
for the Efimovian spectrum for n — >■ +00: 

M)u _ 23/2|sBo|(C - Mmin)(l + C) J,, -(n-l/2)7r/|sBO I 



C2(l + a)i/2 



(108) 



where we have performed the change of variable u 
K{r)r in the integral over r and one has au^-^^/2 = [l 
1/2)2. rpj^g quantity Jd is defined by the integral 



Jcl = -1 + 



C 



du 



1 



Me'' 



fu^ 



VC2 



1/2 



C-w 



(109) 



Whereas as expected [53| this semi-classical result is dis- 
astrously bad for a close to the critical value Uc (it does 
not predict a finite (Zgiob^i at the Efimovian threshold), 

it becomes increasingly accurate for increasing a — ac . 
Comparing to the exact numerical values of FigH^, we 



found that for a 



,m 



-W 



is ai- 



de • > 8, the error on q^i^^^i 
ready less than 10%. In the limit a — )► -|-oo, the semi- 
classical result allows to recover exactly the quantum re- 
sult Eq. (j93p . since Mmin — >■ in that limit and Jd then 
tends to J of Eq. ([M)) . Keeping the first correction linear 
in Uniin hi the semi-classical result we even get the refined 
estimate 



(0 

S'global 



i?* 



2(1 + C)e-^ 



|sBo| 



0(lnQ;/a) 



(110) 

As this amounts to keeping the first a ^/^ term in a 
large a expansion, we can again hope that the Born- 
Oppenheimer approximation (combined to the semi- 
classical one) gives the exact result in Eq. (IllOp ^41]. In 
Fig|4j: we have calculated Jci and hence the semi-classical 
value of ^gioijai numerically to show how it nicely interpo- 
lates between the large-a exact data of Fig|4^ (that still 
strongly depend on I) and the /-independent a -^ -foo 

limit of Ciobai- Taking values of a as large as in FigHJ: 
is straightforward in the semi-classical formula, but it 
would be a numerical challenge for the exact expression 
Eq. (j73p (not to mention real experiments). 



V. NUMERICAL SOLUTION FOR TRIMER 
STATES 

In this section, we proceed with the direct numerical 
solution of the integral equation ([M]) , looking for the al- 
lowed bound state energies E for various values of the 
angular momentum quantum number I. The motivation 
is to look for trimer states that are not predicted (or 
not faithfully predicted) by the analytical results of sec- 
tion IIVI First, in presence of the Efimov effect {I odd, 

a > ac ) the analytical formula ((55)1 is guaranteed to 



be asymptotically exact in the large quantum number 
n — !■ +00 limit, but the numerics can assess its accuracy 
for low values of n, n > 1, and can check whether or 
not n = 1 in (|86p corresponds to the ground state trimer 
for a given angular momentum. Second, it is in princi- 
ple possible that the narrow Feshbach resonance model 
exhibits for 1/a = non-Efimovian trimers, that would 

appear for a mass ratio lower than ac . In single channel 
models, with real interaction potentials, such few- body 
bound states were recently observed numerically [55l . |56[ 
and their emergence was related to few-body resonances 
[56| that one may expect within the zero-range model 
when Eq. ([36]) has a real root s between and 1 [23, [SJ] . 
To be complete, let us mention that, for a Bnite and 
positive value of the scattering length, we expect that 
there exist a finite number of non-Efimovian I = 1 trimer 
states for mass ratios below the critical mass ratio ac~ , 
as shown in [58]. These interesting trimer states should 
have a vanishing energy right at the Feshbach resonance 
(1/a = 0) so that they shall not show up in our numer- 
ical solution and their study is beyond the scope of the 
paper. 



A. Optimized numerical method 

The general method to numerically find the bound 
state spectrum is to approximate the operator appearing 
in the momentum space integral equation by a matrix, 
after discretization and truncation of the momentum fc, 
and to perforin a dichotomic or Newton search of the val- 
ues of the energy E < such that the resulting matrix 
has a zero eigenvalue. Here, we are in the particular case 
of an infinitely narrow Feshbach resonance with an infi- 
nite scattering length, and a much more direct method 
can be used. The kernel in (I34p does not indeed involve 
any interaction length, it exhibits hq as the only momen- 
tum scale, where the wavenumber q is the unknown since 
E = —h^q^/{2ii). We can thus rescale all wavenumbers 
by q, setting k — k/q, K — K/q, and we introduce a 
reduced relative wavenumber 



qrc\{k) 



qrciiqk) 



l + 2a 

(l + a)2 



-,1/2 



(111) 



For convenience, we also write the unknown function 
/«(fc) as 



f(^){k = qk) = ^ 



PH'k) 



kqvciik) ' 



(112) 



where the denominator shall ensure that the resulting 
integral operator is hermitian. After multiplication of 
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34| by k/q,:oiik), we obtain 



qR*r>ik) 



1 



Qrclik) 



P\k) 



du ■ 







Pi[u)kk 



7^ 9rcl(fc)grol(^) J-l I + k"^ + K^ + 



-^kku 

l-\-a 



(113) 



Remarkably, the dimensionless quantity —qR^ is simply 
the solution of an eigenvalue problem for a fixed opera- 
tor. After numerical discretization and truncation, one 
simply has to diagonalize once a real symmetric matrix, 
which is a well mastered numerical problem, and each 
negative eigenvalue of that matrix will provide a numer- 
ical approximation of the quantity —qR* for the bound 
states. We expect the numerical truncation to be accu- 
rate if the maximal value fcmax of k in the numerical grid 
obeys, for each considered negative eigenvalue —qR^: 



^"'max^* ^ Ij 



(114) 



so as to ensure that the effective range term in the two- 
body scattering amplitude is well taken into account. For 
a given numerical diagonalization, the negative eigenval- 
ues —qR* that are larger that — 1/fcmax thus can not be 
trusted. In these estimates, we have dropped for sim- 
plicity a possible dependence of the criteria on the mass 
ratio a l59i1 . 



In presence of the Efimov effect, the condition (|114p 
is quite severe, as gi?* may assume extremely small 
values. The way out is well known, one simply has 
to use a logarithmic scale, with the change of variable 
X = \nk,X — \nK. To keep the hermiticity of the oper- 
ator, we reparametrize the unknown function 



PHk = en 






(115) 



Multiplying (|113p by k^^^ and performing these changes, 
we finally obtain the numerically useful form |60j]: 



-qR,F^^\x) = 



F(')(a;) 



(1 



^2X (.Qg2 jy\ I 



+ OC 



dx 



F^^iX) 



-oo TT [(l + e2^cos2i/)(l-He2^cos2i/)]^/2 



du — 
1 1 



Pi(u) e3(-+^)/2 



g2a; _^ g2X _^ 2ue^+^ siu I 
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r 



where we recall that v = arcsin[a/(l -I- a)], see (|39p . An- 
other interesting feature of the form (I116P , that it shares 
with (|113p . is that the kernel remains bounded at large 
momenta. 



B. Efimovian results 

We now present results obtained from a numerical solu- 
tion of the eigenvalue problem (|116p . for not too large val- 
ues of the mass ratio so that the trimer spectrum does not 
exhibit the hydrogenoid character predicted in subsection 
IIVCI In the numerics, we took for the variable x — \x\k 
a discretization with a constant step dx — 0.09 over an 
interval from Xmin = In (10~^) to a;,„ax = In (10^^). We 
explored the interval of mass ratio a = m/M from the 
small value 1/200 to the large value 200. 

We first explored the case of even angular momenta. 
According to the zero-range theory, numerical obser- 
vation of bound states in that case would reveal non- 
Efimovian trimers. We went up to an angular momen- 
tum I = 12 without finding any bound state: The min- 
imal eigenvalues —qR,, found numerically were positive 
and of the order of 10~^^ or larger. We then explored 
odd angular momenta, looking for non-Efimovian bound 

states for a mass ratio a < ac (and a < 200). We went 
up to / = 13 without finding any. 



At this stage, it remained to explore Efimovian physics. 
The numerical results for the most bound trimer states 
are shown for / = 1 in FiglBl In FiglB^ we directly show 
the obtained values of qR^ as functions of the mass ratio 
a, obviously in log scale for the vertical axis. This is 
useful to estimate which trimer states may be accessed 
in an experiment. An experimental limitation is that the 
scattering length a is not infinite, which will suppress 
the too-weakly bound trimers, that is the trimer states 
with a too small value of q. According to (|28|) we see 
that if l/|a| <C q, the term 1/a is small as compared to 
the term qi-ci{k) for all values of fc, so that the assumption 
1/a = should be a good approximation. Assuming that 
producing in a controlled way a scattering length larger 
than 100 /xm (in absolute value) becomes unrealistic, in 
particular for a narrow Feshbach resonance, we thus take 



,cxp| 



< 100 Aim 



(117) 



and we impose q > 10/|a| « 10/(100^m). Since i?* 
should be much larger than the typical van der Waals 
length (typically of a few nanometers), we take i?* > 
10 nm, so that gi?* > 10~^. Hence the interval of values 
on the vertical axis in FigjS^ [61[ . 

In the more theoretical Fig|6}3, we show the ratios 
of the numerical values of qR^, to the values obtained 
from the approximation (j86p . as functions of the mass 
ratio. The numerical values are numbered as n = 1 
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(ground state trimer state), n = 2,3,..., for q in de- 
scending order (that is in ascending order of the ener- 
gies En — —h^q^/2^). The numerical g„ are then each 
divided by the approximate (/„ of ([55)1 . For a close to 

ttc , for example a < 20, it is found that the ratio of 
numerical to analytical values of the g's is extremely close 
to unity: For the ground trimer ti = 1, the deviation is 
less than 10~^, and for n = 2 and n = 3 it is fluctuat- 
ing between ±10~^ which is probably due to numerical 
errors. Similar results are obtained for larger odd values 
of / (not shown) . 

We have thus numerically obtained the important re- 
sult that the quantum number n = 1 in the analytical for- 
mula (|86p indeed corresponds to the ground trimer state 

(within each subspace of fixed /). For a — ac not too 
large (away from the hydrogenoid regime), the binding 
energy of the ground trimer state thus does not corre- 
spond to the naive expectation « /i^/(/ii?^), it contains 
an extra factor exp(— 27r/|si|) which can be tiny. There is 
numerical evidence that this also holds for three bosons 
resonantly interacting on a narrow Feshbach resonance 
[41.] : the impossibility in the bosonic case to connect the 
spectrum continuously to a S*; = limit (and to study 
the variation of the three-body parameter in that limit) 
however makes the statement more subjective than in the 
present fermionic case. 

To be complete, we have also calculated numerically 
the eigenvectors x -^ F'-'' {x) corresponding to the eigen- 
values —qnR*, in the case I = \ iov a 
cording to Eqs. (|112lll5p one has F^'-^x) 
e^^ cos^ iy)i/2^(0 (qe^). For k/Q < 1, that is e^ cos z/ < 1 
since Q ^ q/ cos v, we found that F^^'^^^ {x) oc e^^/^ as 
deduced from ([54)1 [multiplied by i]. For kR^ cos v ^ 1, 
that is e^cost^ ^ e'^^/'^^l, where si = iSi is the purely 
imaginary Efimov exponent for / = 1, we found that 
F('=i)(a;) oc e'^^/^ as predicted by [42l|. Finally, in the 
crucial intermediate region q/ cosv < k < l/(i?* cosi^), 
which corresponds to the matching interval of the {E < 
0, i?, = 0) and (F = 0, i?« > 0) analytical solutions, see 
Fig 121 we compared the numerics to the analytical result 
deduced from Eq. (l55l) [multiplied by i]: 



14. Ac- 
e3./2(i + 



F 



(1=1)1 



(a;) oc e^' sin[|si|(x — a;o)] 



(118) 



where xq — (arctan |si|)/|si| — ln(2cosi'). From the nu- 
merics, see FiglTl we found that x — > e~^/^F('^^^(a;) 
indeed exhibits half an oscillation of the sinus for the 
ground trimer n ~ 1 [the function remains positive ev- 
erywhere] , and a full oscillation of the sinus for the first 
excited trimer n — 2 [the function changes sign once]. 

These nodal properties were expected from the fact 
noted in 60] that the matrix elements of the kernel in 
(|116l) are strictly negative for I odd, for all x and X, 
whereas the diagonal element is strictly positive, for all 
X. From a standard variational argument, that formu- 
lates Eq. (|116l) in terms of the extremalization of an "en- 
ergy" functional for a fixed norm, and that compares the 
"energy" of a; — )► F*^')(x) to the one of x ^ |F(')(x)|, we 



conclude that the function F^'^ (x) for the ground trimer 
state has a constant sign. As the other modes have to 
be orthogonal to the ground mode, this shows that the 
ground trimer state is not degenerate, and that the ex- 
cited trimer states have a sign-changing function F'^'^(a;). 
We have thus reached a fully consistent picture of the fact 
that n = 1 in Eq. (1861) is indeed the ground trimer state. 



C. Hydrogenoid results 

We now explore numerically the trimer spectrum for 
extremely large values of the mass ratio a. In this limit, 
the low-energy trimers are not expected to be Efimo- 
vian anymore: According to the Born-Oppenheimer ap- 
proach of subsection IIVCI the spectrum for fixed values 
of the quantum numbers n and I becomes hydrogenoid 
for a — )> -|-oo. To test the asymptotic analytical formula 
(|105l) , we plotted in FiglSl the ratio of the wavenumber 

qn [such that E„' — — ?i^[q„ ]^/(2^)] to the asymptotic 
prediction (|104p as a function of l/a^/^, for / = 1 and a 
few values of n (we recall that n — \ labels the ground 
trimer state for fixed X). The prediction Eq. (llOSp then 
corresponds to straight lines, and we indeed observe that 
the numerical results approach these straight lines for 
diverging a. This suggests that (|105|) . obtained in the 
Born-Oppenheimer approximation, is actually asymptot- 
ically exact. As expected, for increasing n, the trimers 
become spatially more extended, so that a larger value 
of a is required to make them hydrogenoid. 

Some considered values of a in FiglSl are extremely 
large, up to 10^. This can not be realized with atoms, 
as such a large mass ratio does not exist in the periodic 
table. Using an optical lattice as suggested in [63 is pos- 
sible if the trimer states have a spatial extension much 
larger than the lattice spacing. This may require here 
unrealistically large values of i?,. A futuristic alterna- 
tive is to replace the fermionic atoms by large and round 
molecules, cooled to their internal (vibrational and rota- 
tional) ground state, that have a vanishing total angular 
momentum in that ground state, may be in the class of 
fullerenes 63| , and that would exhibit a narrow Feshbach 
resonance with the extra atom. 



VI. CONCLUSION 

We have performed a detailed study of the quantum 
three-body problem of two same-spin-state fermions of 
mass m interacting in free space with a distinguishable 
particle of mass M on an infinitely narrow Feshbach 
resonance, with a focus of the three-body bound states 
(trimer states) of that system. The interaction was as- 
sumed to be tuned right on resonance, with an infinite 
s-wave scattering length a, which makes it possible to 
obtain analytical results, since the only length scale left 
in the problem is the so-called Feshbach length i?* [7[. 
The assumption 1/a = also ensures that there is no 



18 



,, 0.1 r 



0< 



0.01^ 



0.001 




100 120 140 160 180 200 

a 




200 



FIG. 6: Numerical study of the Efimovian part of the spectrum: For an angular momentum I — 1, quantities qn (labeled 
in descending order) from n — 1 (ground trimer state) to n = 8, as functions of the mass ratio a = m/M. We recall that 
the trimer energies are then E„ — — /i^g^/2/i. In (a) the numerical values are shown in log scale, restricting the figure to 
experimentally accessible values qRt > 10~^ (see text). In (b) the absolute value of the deviation from unity of the ratio of 
the numerical q„ to the analytical approximation (1861) . The solid (respectively dashed) lines correspond to numerical values of 
qn larger (respectively smaller) than the analytical ones. For both (a) and (b), the vertical axis is in log scale, and the curves 
n = 1 to n. = 8 are from top to bottom. 



1 


' : V 


_^ 1 


,-1- 










w 


; / 


\ 


n=l 


^ II 














H^ 










\ 


/ 








/ 








/ 






\ 


/ 










u 


/ 




\ / 










-1 


'' 1 ■ 


1 


^--'l 1 


1 


■ /^ 


~^ 


A^ 










X 


■ / 


\ n=2 


\ 


^ II 


; / 




\ 


.-^ 


\ ' 




I 


\ '/^"^ 






















f 








/ 






u 




\ / 


- 


-1 




n 


1 



10 



20 



x=ln(k/q^'^) 



30 



0.95 



^^ 0.9 

00 

^ 0.85 

* 

1 ^ 0.8 
0.75, 



^^^^^z 


1 , 1 


^v^^S'-^r^'^"'^-^ 


n=l 


^XV>v~ 


^-^^.y 


; ^>^ 


'S^^ 


n=4 




- 


^V • ^\ . 




N. • 


Ill, 


N,^ • 



0.0025 0.005 0.0075 0.01 0.0125 

1/a 



FIG. 7: For I = 1 and a mass ratio q = 14, ground eigenvector 
(n = 1) and first excited eigenvector (n — 2) of the eigenvalue 
problem Eq. (|116p . corresponding to the ground and first ex- 
cited trimer states. Solid lines: Numerical results. Dashed 
lines: Analytical form (|118[1 . meaningful over the matching 
interval 1/cos;/ < k/q < l/{qR* cos;/) whose meaning is ex- 
plained in FigO and whose borders are indicated by the ver- 
tical dotted lines. The expression of D{k.) in terms of F^^'{x) 
and yi'(k) can be obtained from the text, with x — [n{k/q). 
An overall factor exp(— a;/2) was applied to F^^'{x) for con- 
venience, and the resulting functions are normalized to the 
maximal value of unity. 



two-body bound states. This three-body problem how- 
ever remains rich, richer than e.g. the problem of three 
resonantly interacting bosons on a narrow Feshbach res- 
onance [3, [m, [I4I , because there is a tunable parameter 
left, which is the mass ratio a = m/M of a fermion to 
the extra particle. The existence of this tunable param- 
eter raises the following three fundamental questions on 



FIG. 8: (color online) Numerical study of the hydrogenoid 
part of the spectrum: For an angular momentum 1 = 1, 
wavenumbers g„ for n = 1 (ground trimer state), n = 2, 
n = 3 and n = 4 (from top to bottom), divided by the hy- 
drogenoid asymptotic Born-Oppenheimer prediction (|104[) . as 
functions of 1/q^' . Disks: Numerical results. Straight lines: 
Asymptotic result (|105p including the first deviation of the 
Born-Oppenheimer potential from the Coulomb form at short 
range. The color code is black for n = 1, red for n = 2, green 
for n = 3 and blue for n = 4. 



the trimer states within each sector of fixed total angular 
momentum I. 

First, does this system support trimer states for a mass 
ratio a smaller than the minimal mass ratio ac required 
to activate the Efimov effect |J| ? Such non-Efimovian 
trimer states may indeed emerge from three-body res- 
onances recently discovered numerically for a different 
interaction model [55l |56|, with an Efimov exponent s 
having a rea7 value between and 1 [23, [53 • For the 
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narrow Feshbach resonance, our numerical answer to this 
question is negative. 

Second, how does the Efimov trimer spectrum emerge 
when the mass ratio a is varied across the critical value 
Oc (necessarily for I odd) ? This is an intriguing ques- 
tion, because there is no trimer state for a < a^ (for 
fixed /) and there is an infinite number of trimer states 
for a > Qfc- We found that, for a tending to Uc 
from above, the whole trimer spectrum, including the 
ground trimer state (and not simply the trimer states 
with large quantum number n) forms a geometric se- 
quence. This of course can not be deduced from gen- 
eral zero-range Efimov's theory, which guarantees the 
geometric nature of the spectrum only in the asymp- 
totic region of large quantum number n, and can not say 
anything about e.g. the mo del- dependent ground trimer 
state. We have also shown that the global energy scale 
^global in that spectrum has a finite and non-zero limit 
for a ^ ac- This simply means that, at the Efimovian 
threshold, the energy of the ground state trimer vanishes 
as -Bgiobaiexp(— 27r/|s|), the energy of the first excited 
trimer vanishes as i^giobai exp(— 47r/|s|), and so on, where 
the modulus of the purely imaginary Efimov exponent s 
vanishes as the square root of a — ttc. This constitutes 
a complete picture of the emergence of the Efimovian 
trimer states when the mass ratio is varied across (Xc- 
The dependence of -Bgiobai on the mass ratio a and on the 
angular momentum / was further studied analytically, by 
an efficient and to our knowledge original expression of 
£^giobai in terms of an asymptotic series, see Eqs. (|87l89p . 

Third, what is the nature of the most bound trimer 
states, for example the ground state, when the mass ra- 
tio a becomes significantly larger than the critical value 
ttc? As already mentioned, Efimov's theory can not an- 
swer this model-dependent question. For the narrow Fes- 
hbach resonance, we found that the low energy part of 
the spectrum becomes asymptotically equivalent to an 
hydrogenoid spectrum for a diverging mass ratio, that 
is scaling as -Eo/(" + O^j where n > 1, Ms odd, and 
-Eo = — m?i^/(4Mi?*)^, and we calculated the first devia- 
tion from this hydrogenoid spectrum for a finite a. As the 
hydrogenoid nature asymptotically takes over the Efimo- 
vian nature of the trimer spectrum for a — > -|-oo (except 
in a vicinity of the i? = accumulation point where the 
spectrum remains geometric), in a continuous way when 
the mass ratio is varied, this constitutes a crossover from 
an Efimovian to a hydrogenoid spectrum. 

We have also discussed to which extent all these pre- 
dictions for the 2-1-1 fermionic problem on a narrow Fes- 
hbach resonance may be addressed experimentally with 
cold atoms, the variation of the mass ratio being obtained 
by combining the discrete tuning provided by the choice 
of appropriate species for the fcrmions and the extra par- 
ticle, with an additional continuous fine tuning of the ef- 
fective mass with an optical lattice j62|. Reaching the 
large mass ratios required to observe the hydrogenoid 
part of the spectrum is challenging, except for the fu- 
turistic alternative of replacing the fermionic atoms with 



massive and round molecules. On the contrary, the Efi- 
movian effect in our system may be reachable experi- 
mentally for angular momentum Z = 1, where the critical 
mass ratio is only Oc — 13.607, provided that the energy 
ratio exp(27r/|s|) is not too large, that is the mass ratio 
is far enough from the critical value ac where \s\ — 0. 
From the bound (|117p on the achievable s-wave scatter- 
ing length a with a magnetic Feshbach resonance, the 
ground state trimer is directly observable for a > 15, 
whereas directly observing at least two Efimovian trimer 
states, to check the geometric nature of the spectrum, 
requires a > 20. On a Feshbach resonance as narrow 
as the one of ^Li with ^°K, however, the bound (|117p 
is probably too optimistic, if one does not implement a 
magnetic field stabilization of metrologic quality [6l| . 
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Appendix A: Zero-range solution in position space 

In the case i?* = 0, our three-body problem reduces 
to the zero-range infinite scattering length problem that 
Efimov solved in real space [J] with the ansatz for the 
three-body wavefunction: 



■0(ri,r2,r3) = 



FjR) 



(/7(a2i 



sin(2a3i) 



sin(2a2i 



vmi 



ra 



ar2 -I- ri 



V2 



ara + ri 



1 



a 



(Al) 



This contains two a priori unknown functions of sin- 
gle variables, F and (p. The various variables in this 
ansatz naturally appear in appropriately normalized Ja- 
cobi coordinates. Using e.g the convention of Appendix 
3 in .29j (except for a permutation of the indices), tak- 
ing as in Eq. (fT3|) particle 1 as the extra particle of 
mass nil — M, and 2 and 3 as the fermions of mass 
TO2 — r7i3 = TO, we define the unnormalized Jacobi coor- 
dinates yi = r^ — {mr2+ Mri) / {m+ M) and y2 = ri — r2, 
the generalized reduced masses /ij"^ = m~^ + {m + AI)^^ 
and ^2^ = M~^ + m~^, the hyperradius R > such 
that mR^ = X]i=i "^i(i"i ^ CI)^ where to is an arbitrary 
unit of mass and C is the center of mass position of the 
three particles. Then the conveniently normalized Jacobi 
coordinates are u^ = {pi/inY^'^yi with i = 1,2, so that 
R^ = u\ + u\. This also puts Schrodinger's equation in 
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a reduced form, 



3 

E 



2mj 



Ar. = 



-Ac 



E-^^"= 



2(2m + A/) 

7^ (A2) 

The angles 0:21 and 0:31 belong to the interval [0, 7r/2]. 
The clever choice is tana2i — U2/U1 and the similar for- 
mula for a^i obtained by exchanging the role of particles 
2 and 3. Geometrically, this amounts to introducing po- 
lar coordinates {R,a2i) in the plane {ui^U2)- Explicitly 
this gives 



tan a2i 
tan 031 



ri — r2 cosi^ 



|r3 - {av2 +ri)/(l + a) 

ri — r^\ cost/ 
\V2 - {aTz +ri)/(l + a) 



(A3) 

(A4) 



where cos v is given by (|40l) . Since the second (Faddeev) 
component in (jAll) is deduced from the first one by a 
minus sign and the exchange of particles 2 and 3, which 
ensures the fermionic antisymmetry of i/", it suffices to 
calculate the action of A^ -I- Au2 on the first compo- 
nent. This is quite simple in spherical coordinates, since 
R and 021 depend only on the moduli ui and U2, and the 
factor involving the spherical harmonic function is simply 
y;"'(ui). Thanks to the factor i?^ sin(2a2i) = 2uiU2 in 
the denominator, one is left with the operator 9^^ H- 9^^ , 
that is the Laplacian in the plane (ui,U2), which has a 
simple expression in the polar coordinates (i?, a2i). One 
finally finds that the ansatz (jAip separated in hyper- 
spherical coordinates solves Schrodinger's equation if ip 
solves the eigenvalue problem 

- ^"{a2i) + -^ — -f{a2i) = sV(a2i). (A5) 

cos^ a2i 

s is a priori unknown, but the general theory of e.g. sec- 
tion 3.3 in 29] guarantees that it coincides with s defined 
in Fourier space by (|35|). As a consequence, the R de- 
pendence of the ansatz also separates and the hyperradial 
function F{R) is found to solve a Schrodinger equation 
for a fictitious particle in two dimensions experiencing an 
effective 1/R^ potential: 



EF{R) = - — 



2fh 



F"{R) + 1-F'{R) 



2mR^ 



F{R). 



(A6) 
The Efiniov effect appears for s^ < 0. 

To determine the eigenvalue s^, one needs to specify 
the boundary conditions for Lp. For ui —^ Q, ip in gen- 
eral does not diverge. Since sin(2a2i) vanishes in the 
denominator, this imposes 



^(7r/2) = 0. 



(A7) 



2 approach each other for a fixed position of their cen- 
ter of mass C12, that is V2 = C12 — r/(l -I- a) and 
ri = C12 + ar/(l -f a) with r ^- 0, V should behave 
as A/r -I- 0{r) where A depends here on C12 — r^. In 
other words, there should be no non-zero contribution 
behaving as r°. To calculate the r° contribution from 
the first Faddeev component, one has to expand (f{a2i) 
to first order in 0121 • In the second Faddeev component 
one can directly set r = 0, that is Yi ~ Y2 — C12. From 
the angular representation (|39p for the mass ratio, we 
find arctan[(l -I- 2Q;)^/^/a] = ■j — v and 



-^'(0) - (-1) 



Ml 



^) 



sin(2i^) 



(A8) 



To try to recover the function Ai(s), we heuristically mul- 
tiply this condition by 2cost//(^'(0) to pull out a first 
additive term as in (1441): 



Us) 



(-1)' ^^^ 



^) 



(^'(0) svav 



(A9) 



Note that Eq. (jA5l) is independent of the mass ratio, 
it occurs in the problem of three spin 1/2 (same mass) 
fermions and even in the case of three bosons. The solu- 
tion of (IA5|) obeying (jA7l) was given in these contexts in 
[3 2( 1 in terms of a hypergeometric function. Reusing this 
solution 



"^("21) = 
X 2F1 



cos'^"^ a2i 



/ + 1 



Z+l-s , 3 2 
^ ,«+ 2;cos a2] 



(AlO) 



leads to (H51) that was checked numerically to coincide 
with (j44p . so that the question mark can be removed in 
(jA9p . For completeness, we note that still another form 
of (p(a2i), in terms of a finite sum, was given in [30||. 



Appendix B: Efimovian spectrum global scale at 
threshold 



At fixed angular momentum /, we evaluate the global 
scale (/l-iobai in the limit a —> ai . Since si = iSi, Si > 0, 
vanishes in that limit, we can simply expand the complex 
number Zi up to order Si, neglecting terms 0{Sf). The 
most difficult part is the bit iSiCi{—Si), that we rewrite 

using (USD, r{2iSi) = r(i + 2^^0/(2*^/): 



iSiCi{~Si) 



r(i + 2z^or(i - ivQ - iSi) 
2r{i)r{-ivo + iSi) 

Ti-iUn + iSi) T{l-iVn-iSi 



riGN 



^T{1 - iun- iSi) T{-iVn + iSi) 



(Bl) 



For U2 — )■ 0, on the contrary, ip diverges. More pre- 
cisely, the Bethe-Peierls contact conditions for an infi- 
nite scattering length impose that, when particles 1 and 



We first need to figure out to which order the roots 
Un depend on uq = Si (remember that the t;„'s do not 
depend on the mass ratio). To zeroth order, right on 
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the critical mass ratio, m„ = u'^ by definition. Making 
apparent the dependence of the function A; with a, one 
has 

Ai(m„;a) = 0. (B2) 

It is clear on Eq. P5|) that A; is a regular function of a, 
so that A;(s; a) around ac varies to first order in a — ac 



(B3) 



and (jB2p simplifies to 

A,(m„;a«)=0(a-a«). 



We also note that A/(s;a) is an even function of s. For 
n = 0, uo = Si and Uq = 0; since is a double root of 

s -^ Ai{s;ac), one has to expand (IB3|) to second order 
in u„ — w^ to get the leading contribution, and Sf varies 

linearly in a — ac as expected. For n > 0, u^^ ^ 
and m^ is a simple root of s — >■ A;(s;ac ), one has to 
expand (|B3|) to first order in u„ — w^ to get the leading 
contribution, so that 



-n,(0^ = 



= OiSf) forn>0. (B4) 



0(a 



To first order in Si, we can thus replace the w„'s in (jBip 
by their values u^ at threshold (for n > 0). 

The last step is to expand the F functions in (JB1|) to 
first order in Si. For any real positive number a; > 0, one 
has 



Tjx + iSi) 

r{x) 



1 + iSiiJj{x) + 0{Sf) 



(B5) 



where ip{z) — F'(z)/r(z) is the digamma function. 
One has in particular V'(l) = ^7i where 7 = 
0.577 215 664 9 . . . is Euler's constant, see e.g. relation 
8.362(1) in [31|. We recall that, by definition, — m„ 
for n > and —ivn for n > are real positive. Fi- 
nally, the product over n in (|87p may be written as 
r(/ + 1 — iSi)/r{l — iSi) and expanded with the same 
technique. All this leads to ([TO]) . 

To reveal the convergence of the series in ([TO]) , it is 
useful to introduce the function $(a;) = ipix) + ip{x + 
1) - 21na;. Then <i>(a::) = 0(l/a;^) for x — )• +00 since 
^{x + 1) = Ina; + 5^ + 0(1/ x^) according to relation 
8.344 in ^. Expressing V- (a;) -|-V'(a;-|-1) as $(a;)-|-21n2;, 
for X = —iu^ and for x — —iVn, and collecting all the 
logarithmic contributions as the logarithm of the product 

n^, V , that one can relate to the curvature in s = 
of the function Ai{s;a) thanks to ([7^ . one obtains 



\si\ s,^o+ 



-ij{l + 1) + 3v^(l) - ij{^tvo) - V(l - ivo) 



In 



V, 



^afA,(0;a«) 



2cosz/ 



+00 



n=l 

(B6) 



As one expects that Un/vn — > 1 for n -^ +00, with u„ — 
i{2n+l+l), from the supposed convergence of the infinite 
product in (17^ . the sum in (|B6I) is convergent. 



In practice, it is found that the sum over n > 1 in 
(|B6p is so rapidly convergent that it gives a very small 
contribution to the result. A useful, easy to evaluate 
approximation is thus: 



lim„ ^Sbai 



(^ + ^) ^-2^il+l)-4>(l+2)+3i^{l) 



i?, COS ly 



X d'^Ai{0;a 



(0^ 



(B7) 



The relative error introduced by this approximation is at 
most ~ 5 X 10^'* for the values of / (from 1 to 11) that 
we have considered. Note that the approximation (jB7P 
is equivalent to neglect all terms with A; > 1 in ([55)1 and 



« 



to take the limit a ^ a 



Appendix C: Alternative representations of the 
function Ci{S) 

To determine 9giobai of Eq. (|57)l , that is the phase 9i of 
the complex number Zi in Eq. (l88l) , in the large I limit for 

fixed a — ac or in the large a limit for / fixed, the infinite 
product form (1731) for the function Ci{S) is inappropri- 
ate, even numerically. We construct here more efficient 
representations, in the spirit having led to Eq. (|B6I) . Re- 
markably, the last form that we construct does not rely 
on the roots and poles of the function A;. 

We split the function Ci{S) of Eq. ([75]) in parts whose 
phase is easy/difficult to evaluate. Since —iu„ is real for 
n > 0, and S is real, using F(z)* — T{z*) we rewrite the 
factors of ([75]) as 



T{-iS 



i) r{-iS - iUn)T{l - iS - 



r(l + iS — iUn) 



|r(l + iS — iun)\ 



(CI) 

A similar rewriting can be performed on the factors in- 
volving —iVn, n > 0, SO that, apart from an infinite prod- 
uct that is real positive and does not contribute to the 
phase of Ci{S), we identify a hard part Di{S) that is 
the product of factors of the form r(z)F(z -t- 1), with 
z = —iUn — iS in the numerator and z = —iVn — iS in 
the denominator. We then set 

*^^^"^^^SnS^ with f{z) = z\nz~z. (C2) 

Here we are in the case Re z > so that the branch cut 
of the logarithm (on the real negative axis) and the poles 
of T{z) are out of reach. Then 



Di{S) = e 



-^S,(S) 



n 



^(— ZM„ — iS) 
, -^{-ivn - iS) 



(C3) 



with 



l^iiS) = 2 J2 ifi^i^n - iS) - f{-ivn - iS)]. (C4) 
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According to the relation 8.344 in [311, *(z) = l + 0(l/z) 
at large [zj so that one expects that the infinite product 
in (jC3| converges more rapidly than the original form 
(|73p . Furthermore, only the imaginary part of S;(S') is 
required, and it can be expressed as the integral 



ImEz(^) 



dS'RiiS') 



with 



Ri{S)^\n 



Ai{iS) S^ + il + l)"" 
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Sf 



(C5) 



(C6) 



Eq. (|C5P holds for 5 = 0, since — m„ and — iw„ are all real 
positive for n > 0. To check that it holds at non-zero S, 
one takes the derivative of (|C5p with respect to S, using 
f'(z) = Inz and f(z)* = f{z*). One then recognizes the 
function h.i{iS) from its Weierstrass representation ([72|l . 
also using vq = i{l + 1) and uq = Si. 

It is possible to go further and to express the phase of 
the infinite product in (jC3|) in terms of derivatives of the 
function Ri{S). Relation 8.344 in [3l| indeed gives Stir- 
ling's representation of In^(z) as an asymptotic series in 
1/z, and one also has the fcth derivative for fc > 1: 



dS^ 



M- 



-iu„ 



^S)] 



-i''{k-l)\ 



(C7) 



{~iun - iS)'' 

and similar relations obtained by taking the complex con- 
jugate or by replacing u„ with w„. Finally Di{S) = 
\Di{S)\eMmiS)] with 



W(5) = - 



ds'Riis') Yl ^-^^^^R?'''\s), 



k>l 



(2fc)! 



^(C8) 

where the B2k are the Bernoulli numbers and Rl 
stands for the {2k — 1)*'^ derivative of the function Ri{S). 
The useful statement is then that 



T{2iSi)r{l + l-iSi)Til 



iSi) 



is real positive 



where we used the fact that Ri (S) is an even function of 
S", and thus ipi{S) an odd function of S. Minor trans- 
formations then lead to Eq. (|89p . In short, these results 
originate from the lemma: For any a; > and S real, 

lm[\nr{x-iS)+lnr{x+l-iS)] = - / dS' \n{x^+S'^) 

Jo 



k>l 



(2fc)! rfS'2'^- 



where here again the series is only asymptotic. 



Appendix D: Eflmovian spectrum global scale at 
large angular momenta 

As we show here, asymptotically exact expressions of 
Jffiobai ^^^ ^ diverging angular momentum / can be ob- 



tained analytically. This requires an asymptotic determi- 
nation of the function Ai{iS). To this end, the most con- 
venient starting point is the expression for A; in Eq. (|46|) . 

Since cos ly tends to zero for a > ac m the large / limit, 
we rewrite this expression using relation 9.131(2) of |3H 
that expresses an hypergeometric function of the variable 
z in terms of hypergeometric functions of the variable 
1 — z. For S real: 



M^S) 



cosiy 



l + (-l)'sin'jy 



r(i±^) 



r(i-H^) 



1 



2cosj^ 



^ ;i + i + is i + i-is 1 2 

X 2F1 ( , ,-;cos ly 

^ , _, I + iS ^ I ~ iS 3 n 
- 2^1 ( 1 + -^—, 1 + ^— , -, cos^ ,y 



(Dl) 



For a fixed value of Si , and thus considering a as a func- 
tion of Si, we expect the asymptotic expansion in the 
large I limit: 






(D2) 



Using (jDip for fixed S, the calculation of the leading coef- 
ficient fli is straightforward, expressing each 2F1 function 
in terms of its defining hypergeometric series, see relation 
(9.100) in [3l|, and taking the large I limit in each term 
of the series. For example, for any natural integer fc: 



(!±i )2 . . . (!±i + A; - 1)2 cos2'= J. 



-,2k 



(i)...(i-f fc-1) fc! 1^+00 {2k)V 



(i±iy,,_{i±l + k-lfcos^'^„ 



-,2k 



— ^ 



(§)...(§ + fc-l) fc! i^+oo(2fc + l)! 



(D3) 



(D4) 



The sum over fc then generates cosh and sinh functions 
of ai. Also, the Gamma functions in (jDip may be ex- 
panded using relation 8.344 in [31|. This lowest order cal- 
culation gives A;(iS';)/cosi^ = 1 —exp(—ai)/ai -1-0(1//). 
Since A;(i5'/) vanishes (to all orders), one obtains ai = 
exp(— ai) so that 



ai=C 



(D5) 



where C was introduced in (|77| in the Born-Oppenheimer 
context. This technique can be pushed in principle to any 
order. We calculated 02 = —C/2 and 03. The fact that 
Si does not contribute to 02 (and contributes to 03 in 
the form of a term Sf) is due to the fact that (jDip is 
an even function of S, and that it is always the ratio S/l 
which appears in the expansion. Turning the expansion 
(|D2p into an expansion for the mass ratio, we obtain (|80p . 
dHI]) and (IS2D. 

This large- Z expansion technique can even be extended 
to the case where Si/l — 0(1) (which includes both the 
previous case of Si fixed and the new case Si /I fixed). To 
leading order, one finds 5*2 ~ (a — ac )02/2 (as in the 
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Born-Oppenheimer approximation) and 



M{tS) 



cosiy 



exp 



-C 






1/2- 



C 






1/2 



(D6) 



This shows that Ai{iS) is a function of S of width oc I 
in the large I hmit, so that the derivatives in ((89)) tends 
to zero in that hmit, and the contribution is dominated 
by the integral over S. Setting xi = Si /I, assumed to 
be bounded as we said, we then obtain an asymptotic 
expression for the global scale of the Efimovian spectrum: 



ln('7SLai^*/2) 



Im [lnr(l + iSi) + lnr(l + 2iSi)] 



One simply needs an asymptotic expansion of Ai{iS) for 
S large of the order of Si — ^ +oo, which implies that 
S and 1/cosi/ both scale as a^^^. For S — > +oo it is 
apparent that the integral over 9 in Eq. (1431) is dominated 
by the contribution of a small interval ending in 9 ^ v, 
since here s — iS. Approximating sin(s0)/sin(s7r/2) ~ 
exp[5'(6'— 7r/2)], we see that the small interval has a width 
scaling as 1/S. We then Taylor-expand P/ (sin 0/ sin j/) 
around 6 ~ ly wp to second order in {9 — v), and we 
perform the integral over 9 extending the lower bound 
of the integral to — c», which generates an expansion in 
powers of 1/5. We can also consistently expand cos:^ 
and sini/ up to second order in tt/2 — v, since 1/S and 
7r/2 — i^ are of the same order. If one sets e = tt/2 — v, 
this gives 



;— >+oo 



, , 1 , , 9x arctanx; 

-31n/--ln(l+a;f)+l- 



xi 



Si 

— Injc- x +0 - 
xi I 

(D7) 



with the function 



T{x) 



1 


\ 1 


xf-x^ 


[ <mr 



(D8) 



In the case where Si has a fixed value, for I — >• +00, one 
has that a;; ^- and one may approximate J- by keeping 
terms up to order xf and x"^ inside the square brackets of 
(|D8|) . so that T{x) = (1 + C)/2 + o(l). This gives (|92l). 
In the case where xi fixed to a non-zero value. Si diverges 
for I — )■ +00 so that the Gamma functions in (|D7I) may 
be Stirling-expanded, leading to 



ln(<7riobai^*/2) 



xi fix 
I 



In 



xi 



+ 00 



(l+x2)l/2 



/ -Hixf- 

Jo XI 



arctan xi 

XI 

x^)Tix)]. (D9) 



Furthermore, if a;; ^ 1, the first two terms in the right- 
hand side of (|D9I) tend to zero, and the integral can be 
shown to approach J + ln(l -|- C), where J is the integral 
([M)) . Remarkably, one then recovers for q„io^^iR* the 
same estimate as in (j93p , which was obtained with a dif- 
ferent limiting procedure {Si — > -|-cx) for I fixed). This is 
may be not surprising, since Si/l ;^ 1 in both cases, this 
is even obvious in a semi-classical picture, see discussion 
below Eq. (jl09|) . where the large a limit is reached for 
aC^/2 > (/ + 1/2)2 (-.^^hich implies Si > I), irrespective 
of the fact that / is large or not. 



HiS) 



COSJ^ 



-eS 



eS 



1 



1(1 + 1) 



+0 



1 
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(El) 



where we used P;(l) = 1 and P/(l) = l{l + l)/2. 

A first application of Eq. (|E1|) is an expansion of Si in 
powers of e. Since Ai{iSi) ~ 0, one finds to leading order 
Sie = C, where C is given by (|77)). Going to next order 
gives a correction of order e^ to eSi. Expressing e as a 
power series in 1/a from cose = a/(l -I- a) gives (|83|) . 

A second application of Eq. (jEip is the derivation of 

the infinite-mass-ratio limit of q„i^^j^i- To leading order, 
A;(z5')/cosj^ is a function of eS, and so is 



P/(5)~ln 



X 



fj2 



x—eS 



(E2) 



This means that the (2fc — 1)*^ derivatives of Ri in 
Eq. (jCSp scale as e^fe-i g^j^j g^^g negligible. The integral 
over 5" in Eq. (jC8[) is a leading contribution that scales 
as Si as revealed by the change of variable x — eS'. Also 
the denominator in (jC9p contributes with a phase factor 
-- exp(2iS'; In 2). Using ^^ one finally obtains ((M)) . 



Appendix F: First correction to the hydrogenoid 
spectrum 

Within the Born-Oppenheimer framework of subsec- 
tion IIVCI for the hydrogenoid part of the trimer spec- 
trum, we apply the first order perturbation theory to the 

1 /9 

1/^23 term of Eq. p03p that we call here SV. In terms 
of the Bohr radius gq = fi^ /{niee^) = AR^/a, the nor- 
malized hydrogenoid wavefunction is |64| 



Appendix E: Efimovian spectrum global scale at 
infinite mass ratio 



^global 



The results (|C9P or (|89| are quite useful to determine 



for 



-cxD for a fixed angular momentum I. 



^i!Hr,s) 



[n + l)aQ 



-"^\ {n-l)\ Y'^ 



2{n + l){n + 2l)\ 



xe 



-r23/[{'n+l}ao] 



2r23 



ii 



[n + ;)ao 



r2; + l 



2r23 



[n + l)aa 



yr (r23) 



(Fl) 



24 



Here L^ is the usual Laguerre polynomial defined with 
the convention of [3l| (and not with the one of |64|). 
After angular integration and the change of variable u = 
2r23/[{n + l)ao] we obtain for the expectation value oi 5V 
in that wavefunction: 



{SV) 



K^a"'^ 



(n-1)! 



MRl [2{n + l)Y'/'^{n + 2l)\ 

x/ duu''+''/''e-^[L^^tl{u)\ . (F2) 

To evaluate this integral, we use the generating function 
technique of [6J]: We define 



Hx,y) 



duuf^+^e-'^^p{u,x)vp{u,y) (F3) 



where eventually we shall set 7 = 1/2 and (3 — 2Z + 1, 
and where the generating function of the Laguerre poly- 
nomials L^ for fixed /3 is given for \z\ < 1 by relation 
8.975(1) in [3l|: 



-f-OO 



f(i{u,z) = ^L^(m)2 



r!=0 



^-uz/{l-z) 

(l-zj/S+i' 



(F4) 



On one hand, since ip/3 is the generating function, 



iix,y)= J2 ^"y" 



(F5) 



so we need the diagonal terms n = m in this series ex- 
pansion. On the other hand, from the explicit form of 



I{x,y) = T{l + (5 + -f) 



[{l-x){l-y)]-' 

(1 _ a;y)l + ,3+7 ^ 



(F6) 



that it remains to expand in a series of x and y using three 
times (1 - Xr = j:t=o ^'r(g - ,^)/[r{q + l)r(-z.)], ,. 
non-integer, to obtain for a non-integer 7: 



+00 



dwM^+^e"" 



^liN 



E 

A;=0 



r(fc-7) 
r(/c + i)r(-7) 



r(/3 + 7 + n-fc) 
r(n - k) 



(F7) 



For J — 1/2, f3 — 21 + 1, expressing the Gamma function 
of integers and half-integers in terms of factorials finally 
gives Eqs. (|105ll06p [65,]. 
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